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Abstract. We obtain the Strichartz inequality 



f j \u{t,z)\ A dg{z)dt<C\\u{Q)t 1/A 

JO J M ■ ' 



for any smooth three-dimensional Riemannian manifold (M, g) which is asymp- 
totically conic at infinity and non-trapping, where u is a solution to the 
Schrodinger equation iut + ^Ajv/m = 0. The exponent H 1 ^ 4 (M) is sharp, 
by scaling considerations. In particular our result covers asymptotically flat 
non-trapping manifolds. Our argument is based on the interaction Morawctz 
inequality introduced by Colliander et al., interpreted here as a positive com- 
mutator inequality for the tensor product U(t, z' , z") := u(t, z')u(t, z") of the 
solution with itself. We also use smoothing estimates for Schrodinger solutions 
including one (proved here) with weight r _1 at infinity and with the gradient 
term involving only one angular derivative. 
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1. Introduction 

The purpose of this paper is to establish a L space-time Strichartz inequality on 
a class of non-Euclidean spaces, namely smooth three-dimensional asymptotically 
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conic Ricmannian manifolds (M, g) which obey a non-trapping condition. We make 
these concepts precise as follows. 

Definition 1.1. A smooth complete noncompact Riemannian manifold M, pa- 
rameterized 1 by the variable z := (z 1 , . . . , z n ), with metric g :— gjk(z)dz^ dz k , is 
asymptotically conic if there exists a fixed compact set Kq C M and an (n — 1)- 
dimcnsional smooth compact Riemannian manifold (dM, h) — parameterized by 
the variable y := (y 1 , . . . ,y n ~ 1 ), with metric h :— hjk(y)dy 3 dy k — such that the 
scattering region or asymptotic region M\Kq can be parameterized as the collar 
neighbourhood 

M\K = (0, e ) x OM := {(x, y) : < x < e , y e dM} 

for some e > 0, with metric of the form 

(1.1) g jk (z)dz J dz = — + — 2 

with the usual summation conventions; here hjk(x,y) is a smooth function on 
[0, eo) x dM such that hjk(0, y) = hjk(y)- We will often use the coordinate r := 1/x 
instead of x in the scattering region. 2 

The metric in co-ordinates (r, y) has the 'asymptotically conic' form 

g = dr 2 +r 2 h jk {^,y)dy j dy k . 

By analogy with Euclidean space, we call r the radial variable, and y the angular 
variable; the latter variable is often also denoted ui or 9 in other texts; here we shall 
reserve Greek letters for 'cotangent' or 'frequency' co-ordinates and Latin letters for 
'spatial' or 'position' co-ordinates. For V G T Z M, we let |^| 9 ( 2 ) denote its length 
with respect to the metric g. 

In particular we may compactify M to M := M U dM by identifying dM with 
{0} x dM in this co-ordinate chart. If hij is independent of x (so that hij(x,y) = 
hij(y)) we say that M is perfectly conic near infinity. Finally, we say that M is 
non-trapping if every geodesic z(s) in M reaches dM as s — > ±oo. 

It will be convenient to select a somewhat artificial, but globally defined radial 
positive weight (z), chosen so that (z) is equal to r in the scattering region (0, eo) x 
dM and is comparable to l/e in the compact interior region Kq, in such a way 
that the function (z) is smooth and obeys the symbol estimates 

|Vi(z>|<^-<z} 1 --' 

for all j > 0. 

Remark 1.2. The most important example of an asympotically conic manifold is 
Euclidean space (M,g) := (R",<5), with K equal to the unit ball K n :— {z : \z\ < 



We use 2 as the spatial co-ordinate on M, in order to reserve the letter x for the scattering 
co-ordinate x := 1/r. While we phrase some of these definitions in general dimension n, we shall 
focus primarily on the three-dimensional case n = 3. We will use the indices j, k to parameterize 
three-dimensional manifolds M or two-dimensional manifolds dM (reserving i for the square root 
of -1), and the indices a, (3 to parameterize the six-dimensional product manifold M X M which 
will play a prominent role in the argument. 

2 Note that dM can in fact be recovered from the Riemannian manifold M by identifying points 
in dM with pencils of geodesies that remain at bounded distance apart as t — » oo. 
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1}, say, with dM equal to the unit sphere S n 1 = {z : \z\ = 1} with its standard 
metric h, and scattering co-ordinates 

for r > ro := 1. This example is in fact perfectly conic near infinity, and is 
also clearly non-trapping. More generally, any compact perturbation of Euclidean 
space will be perfectly conic near infinity, though it may not be non-trapping. 
Any asympotically Euclidean space (R™, g) with decay estimates |V J (<? — $){ z )\ < 
Cj( r ) will also be asymptotically conic, with dM equal to the standard sphere 
S"™ -1 ; it is quite likely that the smoothness (and decay) assumptions we use can be 
weakened substantially. 

Remark 1.3. Note that a non-trapping manifold must in fact be contractible. To 
see this, we observe that the region K := {x > e} is, for e sufficiently small, a 
manifold with the same homotopy type as M which is convex in the sense required 
to apply Theorem 4.2 of This guarantees the existence of a closed geodesic in 
M provided nj(K) ^ for some j > 0. 

Remark 1.4. We note that the class of metrics defined above is precisely the class 
of scattering metrics defined by Melrose [lfij . with the metric written in a normal 
form in the scattering region due to Joshi and Sa-Barreto |llj . 

Let (M, g) be an asymptotically conic manifold, and suppose u is a smooth 
solution to the (time-dependent) Schrodinger equation 

(1.2) iu t + -A M u = 

on R x M, where u(t, z) is a complex field, u(0) £ C°°(M) is a Schwartz function 
(i.e. it and its derivatives vanishes to infinite order at dM) and 

(1.3) A M := ^pd ^ k d k = Tr 9 V 2 , g := det(g jk ) 

is the (negative definite) Laplace-Beltrami operator (with V denoting covariant 
derivatives). It is well known (see e.g. 0) that if it(0) is Schwartz then there is a 
unique global Schwartz solution to l|1.2|) . but our estimates will not depend on any 
Schwartz norms of u. We may rewrite (|1.2|) as 

(1.4) u t = -iHu, 

where H := — ^A is the Hamiltonian; note that H is positive definite and self- 
adjoint with respect to the real-valued inner product 



(1.5) (u,v)m '■— Re / u{z)v(z) dg(z) 

Jm 

where dg(z) := ^fgdz is the usual measure induced by the Riemannian metric 
g. In particular we can generate the usual functional calculus of H, and define 
Sobolev spaces H S (M) for all s £ R as the class of functions u whose H S (M) norm 
IMI,ff s (M) : = ||(1 + H) s / 2 u\\l2( M } is finite, where of course we use (|1.5[) to define 
L 2 (M). 

The purpose of this paper is to establish the following local-in-time L 4 Strichartz 
estimate: 
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Theorem 1.5. For any smooth three-dimensional asymptotically conic non-trapping 
manifold M, and any (Schwartz) solution to (|1.2|l (or 1)1.4)1 ). we have 

(1.6) / / nt^ztf dg(z)dt<C\\u(0)\\ 4 H1/4{M) . 

Jo Jm 

Here and in the sequel constants such as C are allowed to depend on the Rie- 
mannian manifold (M, g), but not on u. In particular once one proves this esti- 
mate for Schwartz initial data u(0), one can obtain the same estimate for general 
H 1 / 4 (M) data by the usual limiting argument. 

We now contrast Theorem ll.5l with existing results. In Euclidean space (M, g) = 
(R",i5) we have the Strichartz estimate (see e.g. an d the references therein) 

(1-7) IML?LJ(RxR") < Cn,q,r,*\\ u ( )\\n-(R.«) 

whenever 

Ti 2 n n 

(1.8) s< - + -<-; 2< 9 <oo; 2 < r < oo. 

2 q r 2 

Here of course 



w||l?l-(RxR") := (j (Y \u(t,z)\ r dz^j dt 



1/9 



In particular the estimate 1)1. 6|l is the case (n, q, r, s) — (3, 4, 4, 1/4), but only locally 
in time: t G [0,1]. In Euclidean space a scaling argument, replacing u(t,z) by 
u(X 2 t, \z) and letting A — ► oo, shows that the estimate (|1.6|l holds globally in time, 
that the exponent H 1 ^ on the right-hand side of l|1.6|l is sharp, and that we can even 
replace the inhomogeneous Sobolev space iJ 1 / 4 (R 3 ) by the homogeneous ij 1 / 4 (R 3 ). 
We believe that the global-in-time estimate also holds in the asymptotically conic 
nontrapping setting, but we have not attempted to prove it here. 

Strichartz estimates are crucial in analyzing the low regularity behavior of non- 
linear Schrodinger equations, see for instance PP, 0], and the remarks in Section 

El 

In Euclidean space there is also a local smoothing estimate 

(1.9) / / \Vu{t,z)\ 2 dxdt< C K \\u{0)\\ 2 Hl/a(R3) 

for any compact set K C R 3 ; such estimates seem to have first appeared in the 
work of Constantin and Saut Sjolin ^2] and Vega [2S|. (Earlier estimates of 
this 'dispersive smoothing' type, for the KdV equation, date back to work of Kato 
In Euclidean space the estimates 1)1. 7JI . (|1.9)l are not directly related — indeed 
they can be proven independently of each other. However, on other manifolds the 
local smoothing estimate seems to be indispensable for proving optimal Strichartz 
estimates. 

For general manifolds (M,g), the local smoothing estimate is in general false, 
even if we localize in time, if there are trapped geodesies; this was shown by Doi 
However, assuming that M is non-trapping, we have the analogue of (|1.9|l in 
arbitrary dimension: 

(1.10) J j \Vu(t,z)\ 2 g[t) dg(z)dt<C K \\u(0)\\ 2 Hl/2(M) . 
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This estimate is due to Craig, Kappeler and Strauss [2| on asymptotically Euclidean 
space (see also JO])- We shall need this estimate as well as several global versions 
of it. Some well-known versions (at least for asymptotically flat manifolds; the 
generalization to asymptotically conic is straightforward) are 

r 1 r \Vu{t, z)\ 2 , , 
(1-11) I J m {z)1+e 9iz) dgdt<C e \\u{0)f H1/2[M] 

as well as the variant 



u(t,z)\ 



2 



(1-12) / / '7 v Vi+J dgdt<C e \\u(0)\\ 2 H~u H M) 

Jo Jm \ z i 

for e > 0. These are proved in Appendix II. Furthermore, if tq > is sufficiently 
large, we have the Morawetz estimate 

r 1 r \fu(t, z)\ 2 , s 
(1.13) / / 7 -f^dg(z)dt<c\\u(o)\\ 2 H1/2{M) 

JO J(z)>r \ Z ) 

where y := xV y = (z) 1 'V y , which is well-defined in the scattering region, and thus 
on the domain of integration if r is sufficiently large. Note that for the purposes 
of this estimate, we could equally well define y as the orthogonal projection of 
V to the hyperplane in T Z M orthogonal to d r , since xV Vi form a basis for this 
hyperplane. (For simplicity in writing our estimates, we will henceforth take J/ to 
be defined globally, and equal to in the compact region (z) < t*o.) 
In this paper, we prove the following generalization of (|1.13(l : 

Lemma 1.6. Let M be an asymptotically conic non-trapping manifold. Let a^ k {z) 
be symbols of order zero on M , i.e. junctions obeying the estimates 

\VTa jk (z)\<C m (z)- m 

for all m > 0. I/u(0) € i? 1/2 (M) and u is the corresponding solution to \1.2\ . th 



en 



a^ k (z)f j u{t,z)V k u(t,z) 
dg(z) 

M r 



dt < C\\u(0)\\ 2 H1/2{M) . 



This improves upon l|1.13f) since only one of the u-derivatives is required to be 
angular. In flat Euclidean space it was proved recently by Sugimoto |22j . In fact, 
we prove (and need) a stronger version of this result where the symbols a? are 
allowed to depend on an additional parameter; see Lemma 14.31 

The epsilon loss in (|1.11|) . (|1.12|) is not removable, even in R™ (although refine- 
ments with logarithmic losses are certainly possible); semi-classically, this can be 
explained by the fact that a particle z(s) moving at unit speed along a geodesic will 
have a finite integral / R (^(s)) _1_e ds, whereas the integral J R (z(s)) _1 ds diverges 
logarithmically. The Morawetz estimate (|1.13|) then corresponds semi-classically to 
the fact that the slightly smaller integral J/ z f s \\ >ro l w ang(s)| 2 (z(s)) -1 ds avoids the 
logarithmic divergence and converges absolutely, where v ang (s) — r-j^y(s) is the 
angular component of the velocity. This reflects the fact that the velocity vector of 
a geodesic becomes purely radial in the asymptotic limit. 

Now consider the Strichartz estimate (f 1 . T|t on general manifolds. Without the 
non-trapping condition, one does not have any local smoothing, and one does not 
expect to obtain the estimate (|1 . 7|1 with the sharp number of derivatives s. However, 
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it is still possible to obtain a (local-in-time) Strichartz estimate with a loss of 
derivatives. Indeed, Burq, Gerard and Tzvetkov showed that the estimate 

(1-14) IMIl ? l s([ o,i]xm) < C\\u(0)\\ HS+i{M) 

holds for arbitrary smooth manifolds (trapping or non-trapping) in general dimen- 
sion. This loss of derivatives, as well as the localization in time, was shown in |3] 
to be sharp in the case of the sphere (which is in some sense maximally trapping) . 

The next progress was by Stafhlani and Tataru |2(J| . who were able to remove 
the 1/q loss of derivatives for metrics on R" which are non-trapping and Euclidean 
outside a compact set (and also only required C 2 regularity of the metric); a key 
tool was the use of to localize to a compact region of space. After the work 

of Stafhlani and Tataru, Burq gave an alternative proof of the same result, and 
also considered metrics which were asymptotically flat rather than flat outside of a 
compact set. In this more general setting Burq was almost able to recover the same 
Strichartz estimate as Stafhlani and Tataru, but with an epsilon loss 3 of derivatives: 

IMU?££([0,1]XM) < C\W(0)\\h'+*(m)- 

The idea was to divide M into dyadic shells R < \z\ < 2R and apply a variant 
of the Stafhlani- Tataru argument on each shell, either rescaling R to become 1, or 
relying on instead of l|1.10fl . As remarked earlier, the presence of the e in 

(|l.llfl is essential, and is related to the epsilon loss of derivatives in Burq's result. 

Our result in Theorem 11.51 is slightly stronger than Burq's in the sense that 
it applies also to asymptotically conic manifolds, and removes the epsilon loss. 
However, it is also weaker than these previous results because it is restricted to the 
specific exponents (n,q,r,s) — (3,4,4, 1/4) and cannot obtain the full range (|1.8|) 
of the estimate l|1.14[) , although one can obtain a subset of this range by applying 
Sobolev embedding to i|1.6fl and also using the energy identity \\u(t)\\ L ^ L 2^- RxM - ) = 
ll M (0)IU 2 (M)- We conjecture however that the full set of Strichartz estimates can 
be recovered with no loss of derivatives for non-trapping asymptotically conic (or 
asymptotically flat) manifolds in arbitrary dimension. 

Unlike the previous results, the proof of Theorem II . 51 does not proceed via con- 
struction of a parametrix or proving any dispersive estimates on the fundamental 
solution. Instead, we use the interaction Morawetz inequality approach introduced 
in 0, [7] in the context of a non- linear Schrodinger equation in Euclidean space. 
We re-interpret this approach in the language of positive commutators applied to 
the tensor product U(t,z',z") :— u(t, z')u(t, z") of the solution u with itself, and 
then modify it to asymptotically conic manifolds. 

There will be a number of technical difficulties in carrying out the above strategy, 
first in understanding the error terms generated by the non-zero curvature, and 
second the problem of avoiding the singularities of the metric function <Im{z' \z") 
once this distance exceeds the radius of injectivity. Unlike the argument in [5], our 
approach cannot afford to localize the radial variable to dyadic blocks R < r < 2R 
as this would create error terms with radial derivatives in them, which cannot be 
controlled by (|l.ll(l without losing an epsilon too many derivatives. In fact we 
can only localize r once, to the region r > tq, though we can localize the angular 



The endpoint q = 2 is also not obtained, for a rather technical reason involving the causality 
cutoff s < t in Duhamel's formula. 
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variable y to a small set 4 . The understanding of the geometry in this regime — in 
particular, controlling the derivatives of the distance function dM({r', y'), (r" ,y")) 
when y',y" are close, and r',r" are both large but not necessarily comparable — 
appears to be absolutely essential for obtaining bounds such as (|1.7|l with no loss of 
derivatives whatsoever. Thus our argument requires a certain amount of geometric 
machinery which we have placed in an appendix (Section 

It is well known that Strichartz estimates can be used in the theory of nonlinear 
Schrodinger equations, and we briefly discuss a modest application of this estimate 
to a class of non-linear problems in Section[S] However as our Strichartz estimate is 
restricted to L in space and time, the application is admittedly somewhat limited. 
It is still an important problem as to whether one can obtain the full range of 
Strichartz estimates (with no derivative loss) on general non-trapping manifolds, 
and whether they can be extended to be global in time; such a result will have 
substantially more applications to non-linear Schrodinger equations. 



In this section we review the basic method of positive commutators in the context 
of the time-dependent Schrodinger equation, and also indicate how it is connected 
to estimates of Morawetz type; we will rely on these methods, together with de- 
tailed analysis of the various terms arising from such methods, to prove Theorem 
11.51 In this discussion M is a general smooth Riemannian manifold with compact- 
ification M. In our applications M will either be a three-dimensional non-trapping 
asymptotically conic manifold, or else the product manifold M x M which is six- 
dimensional and non-trapping, but not asymptotically conic (as we shall see, it has 
a 'corner' at infinity which we will have to blow up to analyze). To emphasize 
that we may be working in the product setting, we will use Greek indices a, (3, 7 to 
parameterize the manifold rather than Latin indices j,k,l. 

Let u be a solution to the Schrodinger equation (|1.2fl (or l|1.4|l K to avoid technical 
issues let us assume that u e C°°{M), that is, u is rapidly decreasing together with 
all derivatives as z — > 00. Given a pseudo-differential 5 operator A on M, we can 
easily verify 6 the Heisenberg equation 



where [H, A] := HA — AH is the usual commutator. Integrating this we obtain the 



Even doing this localization creates some unpleasant error terms in which only one of the 
derivatives is angular. We shall control such terms using Lemma ll.6l or more precisely Lemma l4,3l 
which only requires one angular derivative. 

^In fact, we shall mostly be able to rely on classical differential operators A for our commu- 
tants; pseudo-differential commutants will only be needed to prove the local smoothing estimates 
necessary to control errors arising from the classical commutant. 

6 Note the presence of the real part in if .51 allows us to obtain this even when A is not self- 
adjoint. In practice, we will only deal with operators A whose principal symbols are real, so that 
they are self-adjoint modulo lower order terms. 



2. Positive commutators and Morawetz identities 



d t (Au(t),u(t)} M = (i[H,A]u(t),u(t)) M , 



inequality 



(2.1) 




8 



ANDREW HASSELL, TERENCE TAO, AND JARED WUNSCH 



The strategy of the method of positive commutators is to select a good commutant 
A, chosen so that the commutator i[H, A] is positive (semi-) definite, possibly mod- 
ulo a lower order error term, and then apply l|2.1[l to obtain a non-trivial spacetime 
bound on u. From the pseudo-differential calculus and the sharp Garding inequal- 
ity, one only expects to be able to do this if the symbol a(A) of A increases along 
the bicharacteristic flow of H (i.e. geodesic flow in phase space). At first glance 
this method only seems able to obtain quadratic estimates on u (e.g. the method 
can be used to deduce (|1.1U|I - l|1.13[) ). but we will also be able to obtain quartic 
estimates (and in particular the L 4 bound Ijl.Gp ) by the trick of replacing u by the 
tensor product U(t,z' , z") := u(t, z')u{t 1 z"), and working in the six-dimensional 
product manifold M x M. 

For future reference we compute some commutators of H with multiplier oper- 
ators. 

Lemma 2.1. Let a(x) be a real-valued tempered distribution on M , thought of as a 
multiplier operator (af)(x) :— a(x)f(x) on Schwartz functions. Then we have the 
single commutator identity 



Proof. The proof of l|2.2jl is elementary. To prove l|2. 3|1 . observe that both sides 
of H2.3JI are second-order self-adjoint operators (with respect to (,)m) with real 
coefficients. Also, by applying both sides to the constant function f we see that 
the constant terms of (|2.3() match. Thus it will suffice to show that the principal 
symbols of both sides match, since this would mean that the difference of the 
left-hand side and right-hand side is a self-adjoint first-order operator with real 
coefficients and no constant term, which is necessarily zero. 

To compute the principal symbols at a point p, we pass to Riemannian normal 
coordinates about p. Then up to second order at p, g a p = 5 a p, and inserting l|2.2[l . 
we see that the LHS of (I2.3|) becomes 



modulo lower order terms. That the symbol agrees with minus the Hessian follows, 
as Christoffcl symbols in (12.4(1 vanish at p. 



Remark 2.2. Because i[H,a] is self-adjoint, we know that J (i[H, a]u,u) dg(z) is 
real for any test function u. Since J (i(Ha)u, u) dg{z) is clearly imaginary, we thus 
see from (|2.2|) that 



(2.2) i[H, a] = i((Va), V) 9 + i(Ha) = i(V a a)V a + i(Ha) 

and the double commutator identity 




(d a df3(a))d a df3, 



□ 



(2.5) 




More generally, if b is another real function, then 



(2.6) 
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Remark 2.3. The Hessian Hess(a) Q ^ measures the convexity of the function a(z) 
with respect to the geodesic flow. Indeed, if z(s) is the curve of a unit-speed geodesic 
in M, and v a (s) := -^z a (s) is the velocity vector, then one has 

(2.7) ^ a ( z ( s )) = Hess(a) Q/3 (z( S )K( S K( S ). 

This can be seen either from <|2.4I) and the geodesic flow equation, or alternatively 
from (|2.3I) and the symbol calculus. 

Example 2.4 (Euclidean one-particle Morawetz inequality; ^3)- Consider flat 
three-space M := R 3 . Let a(z) := \z\, the distance function to the origin. A simple 
computation shows that 

., x z^ z k 

Kess(aY (z) = djdka(z) = — — — — , so that 

\ z \ \ z \ 

Kess(ay k {z)diu(t,z)d k u{t,z) = i^feL 

\ z \ 

where J/ denotes the angular gradient as before. Also, 

-H 2 a = --AAa = 2ttS, 
4 

where 5 is the Dirac mass at the origin in R 3 . Thus by l|2.3l) . I|2.5|l we have 
(-[H,[H,a]]u,u) M = 27r| M (t,0)| 2 + / Iffgjg^L dz 



and 



R :! 



z 

(i[H, a]u, u)m = Ini / u(t, z)-r—r • Vu(t, z) dz. 



\z\ 

Since V maps Hi(R?) to i?-2(R 3 ), and A is a bounded multiplier 7 on H I (R 3 ) 
we have 

< C||u(t)||^ 1/2 ^ R 3^; 



u(t, z)- — r • Vuft, z) dz 

\A 



since the ij 1 / 2 (R 3 ) norm is conserved, setting A — i[H,a] in i|2.1|) we thus obtain 
the Morawetz inequality 



1 r 1 r \V7„.(+ „M2 

2 



u(t,0)\ 2 dt+l I i' u <t> z W dzdt<C\\u(Q)\\l 1/2(R3y 

JO JR 3 1^1 v ' 

In fact, in Euclidean space this particular estimate extends from the time interval 
[0, 1] to all of R by the same argument. Note in particular that we have proven 
(|1.13fl for Euclidean spaces. A similar argument can be used to derive (|1.13(l for 
general non-trapping manifolds, as well as Hl.llf) . if one already knows the dis- 
persive smoothing estimate (|1.10|) : however the estimate H1.10[) cannot in general 
be proven without introducing commutants which are pseudo-differential operators 
rather than classical differential operators. 



7 This can be seen by first verifying boundedness on L 2 (R 3 ) and iif 1 (R 3 ) (using Hardy's 
inequality) and then using interpolation and duality. 
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Example 2.5 (Euclidean two-particle (interaction) Morawetz inequality; |Zj)- 
We now consider a variant of the Morawetz inequality. We let M be the six- 
dimensional space M := R 3 x R 3 , which we parameterize as z = (z',z") where 
z',z" € R 3 and z £ R 3 x R 3 ; observe that the Hamiltonian H := — ^Ajj.3 X r3 on 
this manifold splits as H = H z / + H z " where H z / := — ^A z ' and H z » := — |A Z ". 

Let a : R 3 x R 3 — > R be the Euclidean distance function a(z' , z") — \z' — z"\. 
The function a{z' ', z") is clearly (weakly) convex under geodesic flow of z' and z" 
(as can be seen for instance using center of mass co-ordinates), thus the tensor 
Hess(a) Q ' 3 is positive semi-definite. In particular we have 

(2.8) Ress(a) al3 (z)d a U(t, z', z")^U{t, z', z") > 0. 
One can also compute 

(2.9) -H 2 a = 8it5(z' - z"), 

thus —H 2 a is a Dirac mass on the diagonal {(z' , z") £ R 3 : z' = z"} of R 3 x R 3 . The 
same reasoning as in the previous example yields the product Morawetz inequality 

\U(t,z,z)\ 2 dzdt < C\\U(Q)\\ 2 H1/2(R3xR3y 

This inequality holds for all solutions U(t, z', z") to the Schrodinger equation on 
R 3 x R 3 . If we now specialize to the tensor product U(t, z', z") := u(t, z')u(t, z") 
of an R 3 solution u(t, z) of the Schrodinger equation 8 , we obtain 

(2.10) / f lu^z)^ dzdt<C\\u(0)\\ 2 L2 ^3)^(0)11^ R3 ; 

JO JR 3 v ; 

after a Littlewood-Paley decomposition (see Proposition ^ . II below) . this will imply 
in particular the Strichartz estimate 1)1. Again in the Euclidean space setting it 
is possible to extend the time interval from [0, 1] to all of R. 

We have just proven Theorem 11.51 for Euclidean space R 3 . It is thus tempting 
to adapt the above proof to more general manifolds. The most obvious thing to try 
is to repeat the above argument, with the manifold R 3 x R 3 replaced of course by 
M x M (with the product metric g®g), and with the function \z' — z"\ replaced by 
the distance function duiz' ' , z"). This works well locally, when z' and z" are close 
together; the main issue is to verify that djvr(z', z") remains essentially geodesically 
convex, and that the error terms can be treated by the one-particle estimates (ll.llf) 
- (|1.13|) . Difficulties arise however in the scattering region, requiring more careful 
control on the geometry of the distance function on asymptotically conic manifolds 
over long distances, as well as requiring a refinement of Lemma |l.6l but we shall 
address these issues later. 

In Euclidean space one can also develop /c-particle versions of the Morawetz 
inequality for any k > 2 by working in center-of-mass coordinates, with the function 
a being the standard deviation of the k positions Z\ , . . . , z^. This leads to a useful 



°Note here we are using a very special property of the Schrodinger with itself, namely that 
tensor products of solutions to the Schrodinger equation remain solutions to a Schrodinger equation 
(but now on the product manifold of twice the dimension) . This property is not shared by other 
equations such as the wave or Klein-Gordon equation, and it remains a challenge to extend this 
type of argument to those equations. 
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monotonicity formula in d dimensions when d(k — 1) > 3 (for instance the 4- 
particle Morawetz inequality gives an L\ z estimate in one dimension). In this 
paper, however, we will confine ourselves to the two-particle problem. 

3. Littlewood-Paley reduction 

In the previous section we discussed the method of positive commutators, and 
showed how it can be used to prove Lf x type estimates such as (|1.6|l . It is then 
tempting to apply the same argument immediately to general manifolds, using the 
distance function du(z' ', z") as a substitute for \z' — z"\ as mentioned above. Before 
we do so, however, it is convenient to make a Littlewood-Paley reduction which 
allows us to restrict the solution u to a single dyadic frequency range [2 2k ~ 2 , 2 2fc+2 ] 
in the spectrum oi H. This reduction is standard (see e.g. |19| . Lemma 5.1 and the 
ensuing discussion). 

Proposition 3.1. Suppose we have 

(3.1) [ f \u(z,t)\ 4 dg(z)dt<C2 k \\u(0)\\ 4 L2(M) , 

Jo Jm 

for all solutions u(z, t) of <|1.4|l with u(z, 0) in the range of the spectral projection 
X[2 2fc - 2 ,2 2fc + 2 ] {H)> an d all k > 1. Then Theorem \1.5\ holds. 

Proof. The solution to (|1.4J) is given by u(t) = e~ ttH u(0). Furthermore we have a 
Littlewood-Paley decomposition 



k=l 

where </>o is a compactly supported bump function, and ipk is a bump function 
adapted to the interval [2 2fe ~ 2 , 2 2k+2 ]. This induces a decomposition 

oo 

u(0) = u(0) (0) + u(0) (fe) , u(0) (0) = MH)u(0), u(0) (fe) = ^ k {H)u{0). 

k=l 

Informally, ipk(H) is the spectral projection to functions of frequency comparable 
to 2 fe . 

The contribution of the low- frequency term u(O)^ -' can be handled very easily 
using the Sobolev embedding ||u|| L 4 < C||Vw|| £ 2 in four dimensions: 

\\e- im u(0)W\\ LUlQA]xM) <C\\V z>t {e- itH (f>a(H)u(0))\\ L 2 {[0 , 1]y< M) 
< C\\(l + H)e~ UH MH)u(0)\\ L 2 {[OA]xM) < C\\u(0)\\ L 2 {M) < \\u(0)\\ 
So consider the high-frequency terms 

oo 

(3-2) llE e ^(°) (fc) ll^([o,i]xAf). 

fe=l 

We use the Littlewood-Paley estimate 

\\Y,MH)f\\ LH M)<C\\(Y,\MH)f\ 2 ) 1/2 \\ LH M) 

k k 

(see e.g. or chapter VI, section 7.14; the point is that the heat kernels e~ tH 
or resolvents (H — A) -1 enjoy good decay and regularity properties, and thus the 
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Littlcwood-Paley square function operator is of Calderon-Zygmund type), to bound 
by 

OO 

fc=l 

we can move the I 2 sum outside of the L 4 norm to bound this by 

OO j/£ 

c(Eii e ^(°) (fe) iiit,([o,i]xM)) ■ 

fc=l 

The estimate <|3.1ll implies that <|3.2ll is bounded by 

c(]T2*/ 2 |Ko)«||i 2(M) ) . 

k=l 

From the spectral localization of u(0)( fe ), we can estimate this by 

oo j, 2 

c(x;ii(i+jj) i/8 «(o)W|ii, (i0 ) . 
fe=i 

Since ^(O)^- 1 = ipk(H)u, we see from spectral calculus that this expression is 
0(||(1 + J H-) 1 /Vo)|| i2(M) ) =°(\H°)\\w/*), as desired. ^ □ 

Henceforth we fix k > and write u — thus u is implicitly localized to the 
portion of the spectrum of H lying in the interval [2 2k ~ 2 , 2 +2 ]. Our task is now 
to prove iJjOJ. Observe from this spectral localization (and the unitarity of the 
operators e~ ttH ) that we have the bounds 

(3.3) c2 fcs |K0)|| i2(M) < \\u{t)\\ H , {M) < C s 2 ks \\u{0)\\ LHM) 
for all s e R and t G [0, 1]. In particular, we have 

(3.4) \\Vu(t)\\ LHM ) < C\\u{t)\\ HHM) < C2 k \\u(0)\\ L 2 (M) 

(the first inequality being easily verified by an integration by parts). These in- 
equalities give us a significant amount of freedom to estimate various factors of 
u in various Sobolev spaces, which will be very convenient in the argument that 
follows. Indeed, because of these estimates, the location of the derivatives in our 
expressions to be estimated will no longer be particularly relevant; what matters 9 
is the number of such derivatives (which range between and 2, with the terms 
with two derivatives of course being the most difficult) as well as their type (angular 
derivatives will, in general, be easier to handle than radial derivatives, thanks for 
instance to I|1.13|l 1. 

4. Overview of proof of Theorem 11.51 

Now that we have localized u to a single frequency range, we apply the method 
of positive commutators discussed in Section [3 We will in fact use this method 
twice; once for the near region Kq, and once the scattering region M\K$, using 
different commutants for the two regions. 



^Of course, weights such as l/(z) and !/{«') will also be important to keep track of. 
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Let < £ < 1 be a small constant to be chosen later. We use Z to denote the 
quantity 

(4.1) Z:= £ - 2 2 fc H0)|| 4 L2(M) +e / / \u(z, t)\ 4 dg(z) dt. 

JO J M 

Clearly, to prove l|H.l|l . it suffices to show that 

(4.2) [ [ \u{z,t)\ 4 dg{z)dt = 0{Z), 

JO J M 

since for e small enough we can absorb the second term in 14.1f> into the left-hand 
side. 

It remains to prove 1)4. 2[1 . We first exploit the compactness of M to work on 
smaller regions. 

Proposition 4.1 (Local L A bound). For every z £ M we have 
[ [ \u{z,t)\ i dg{z)dt<C Z0 , v Z 

JO J B(z ,v) 

if rj is sufficiently small (depending on zq), where B(za,r]) is the ball of radius r\ 
centered at Zq. 

Proposition 4.2 (Scattering L 4 bound). We have 

[ [ \u(z,t)\ 4 dg(z)dt <C ro Z 

JO J(z)>r 

if ro is sufficiently large. 

From these two propositions and the compactness of M we easily obtain 1)4.2(1 
as desired. 

It remains to prove the propositions. Proposition ^. H is the easier of the two and 
is proven in Sectional by applying the positive commutator method on the product 
manifold M x M (with Hamiltonian HmxM = Hz' + H z ») to the product solution 
U(t,z',z") u(t, z')u(t, z"), and with the commutant A noar = A ncar zo „ defined 

by ' 

(4.3) A ncar := ip{z', z")i[H Mx M, d M {z', z'% 

where ip(z' , z") is a smooth non-negative cutoff to the region z', z" S B(zo, rj) which 
equals 1 on the region z', z" 6 B(zq, rj/2). Note that while the first-order operator 
^4near is not quite self-adjoint, it does have real principal symbol and its commutator 
will turn out to be positive definite (modulo lower order terms). The commutator 
[Hm x m > ^4near] turns out in fact to be quite tractable, as most of the error terms 
generated can be treated by the standard local smoothing estimates p.lljl . 1)1.12)1 . 

Proposition 14.21 is proven in a similar manner but involves more technical dif- 
ficulties (in particular, one needs much more precise control on the derivatives of 
the metric function dM near infinity) . The proof applies the positive commutator 
method onMxJlf to the same product function U discussed before, but with the 
commutant A asympt = A asymp t, ro now given by 

(4.4) ^asympt = x(z' , z")<ip(y' , y")i[H M xm , d M (z' ', z")} 

where x( z ' \ z ") is a smooth non-negative cutoff to the product scattering region 
{(z'), (z") > r /2} which equals 1 on the slightly smaller region {(z'), (z") > r }, 
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and ip is a smooth non-negative cutoff 10 to the region {ddM(y' ,y") < 2^} which 
equals one on {ddM{y', y") < ??}, where ij > is a small number (less than a quarter 
the injectivity radius of dM) to be chosen later. Note that while x4> localizes y' 
and y" to be close to each other, it does not localized r' and r" to be similarly close 
(although we have made them both large); it turns out that we cannot afford to 
localize further as this will generate error terms that are too difficult to estimate. 

In the asymptotic regime, when controlling [H, A asympt ], terms arising from 
derivatives hitting cutoff functions are harder to control, mainly because there 
are many terms which decay only like l/(z') or l/(z"), and the estimates <|1.II[1 . 
(|f .I2f> instead require decay of l/(z') 1+£ , l/(z") 1+£ . One can use Ijf.f 3fl to con- 
trol those terms which have two angular derivatives. Unfortunately there are also 
mixed terms involving exactly one angular derivative; these terms come from one 
derivative hitting the angular cutoff ip and one hitting the distance function d,M- 
A typical such term has the form 



a ^^yxMQv^(M0 Kti 0|a dgl dg „ ^ 

">r /2 r 

where a y (z',;z") are symbols of order in z', uniformly in z" (see Proposition |^| 
below), and J/ i represents a component of the angular derivative ^V y . If both 
derivatives were angular then such a term could be controlled by the single particle 
Morawetz estimate (|I.13|) . however this estimate is insufficient when one of the 
derivatives has a radial component. Fortunately these mixed terms can still be 
handled by the following new endpoint smoothing estimate: 

Lemma 4.3. Let W be an arbitrary index set (possibly infinite). For each w £ W , 
let a(^(z) be a tensor-valued symbol (j,k range over the co-ordinate indices of dM 
and M respectively) supported on the scattering region {(z) > ro} for sufficiently 
large ro which are uniformly symbols of order in the sense that 

\V™ai k (z)\ < C m (z)- m 

for all 11 m > 0, where the constants C m do not depend on w. Then we have 



f L f a^lz)JMt,zW k u(t,z) 
(4.5) / sup / > r > \±l dg{z)d t<c\\um\ 2 H ^ {M y 

JO w£W Jm r 
We also have the variant 

l.(zW,u(t. z)u(t. z) 

dt < C||u(0)||| a(M) 



(4.6) / sup ( 

JO wGW Ja 



<(z)fMt,z)u(t,z) 

dg\z) 

im r 



where a° w obeys similar bounds but without the additional index k. 

This lemma appears to be new even when If is a singleton set, and may be 
of independent interest. Morally speaking, it corresponds semi-classically to the 
fact that J^ >rn \ v img( s )\/(z(s)) ds converges for all geodesies z(s) of unit speed 
(cf. the discussion after Lemma Tl. 611 : however the result is still rather subtle, since 
the inequalities fail if the absolute value signs are brought inside the z-integral; see 
Proposition 18. II The proof is somewhat technical, requiring use of the scattering 



10 It may seem more intuitive to localize y' and y" separately rather than to localize 
ddM (j/'i ?/'): as ^ turns out we cannot afford to let this cutoff be non-constant on the di- 
agonal z' = z" and so we must proceed in this fashion. 

t^In fact, only finitely many m are necessary. 
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pseudo-differential calculus, and is deferred to an appendix ( Sect ion ITU)l . We remark 
in the case where M is asymptotically flat rather than asymptotically conic, one 
can avoid the need for this Lemma by choosing a more sophisticated commutant 
(see the remarks in Section [&J. 

Note that in our applications of Lemma f4. 31 the symbols a(^(z) and a 3 w (z) will 
often be derivatives of cutoff functions and of the distance function on M x M, with 
w corresponding to one variable and z to the other. 

It remains to prove Proposition 14.11 and Proposition 14.21 This will be done in 
the next few sections. We first however record a number of quantities which are 
bounded by Z. 

Lemma 4.4. We have the estimates 

(4.7) / / (\Wu(t,z')\ 2 \u(t,z")\ 2 + \u(t,z')\ 2 \Wu(t,z")\ 2 ) 

JO JMxM 

X{ W + j^) d 9(z')dg(z")dt = 0(Z) 

(4.8) f 1 f \fu(t,z')\ 2 \u(t,z")\ 2 ± J dg(z , )dg(z")dt = 0(Z) 

Jo J{z'),{z")>r r 

(4.9) [ f \u(t,z')\ 2 \u(t,z")\ 2 (l+ ) )dg{z')dg{z")dt = 0{Z). 
Jo Jmxm <J-m\z , z ) 

(4.10) 



aP k (z'. z")fMt, z')V k u(t, z') 

— } * 3 V , ; - - \u(t,z")\ 2 dg(z')dg(z") dt = O(Z) 



J M x M 



(4.11) 



1 f a jk (z',z")fMt,z'Mt,z') „ 

/ — — 1 / u(t, z")V k u(t, z")dg{z')dg{z") dt = O(Z). 

JMxM r 

whenever a^ k (z', z") are symbols of order in z 1 uniformly in z" that are supported 
on the scattering region {(z') > r } for sufficiently large. 

Proof. We begin with (|4.7(l . Consider for instance the term 

Vu(z')\ 2 \u(z")\ 2 ^dg{z')dg{z")dt. 

JMxM \ z / 

Applying l|l.ll|l in the z' variable and 13.3|l in the z" variable we can bound this ex- 
pression by O(\\u{0)\\ 2 H1/2(M) \\u(0)\\l HM) ), which by ^ again is O(2 fc ||u(0)||| 2(M) ) 
F = O(Z). Now consider the term 
™i 



Hz'^lVuiz'^j^dgiz^dgiz'^dt. 



I^XIV^'Ft 1 

/0 JMxM V 

Now we use (|1.12() in the z 1 variable and (|3.4|l in the z" variable to bound this by 
O(||«(0)||^_ 1/2(M) 2 2fc || W (0)||2 2(M) ), which by O is again O(2 fe || U (0)||| 2 ) - O(Z). 

The estimate (14 .8|) is proven similarly to l|4.7|l . but with l|1.13|l used instead of 

Now consider (I4.9|) . The integral on the region du(z' ', z") < e is clearly bounded 
by 0(e _2 ||/|ll2 (M) ) = O(Z) thanks to (E3J). When d M (z', z") < e we observe that 
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the kernel (cIm(z' , z"))~ 2 has marginal integrals 



sup/ (d M (z',z"))- 2 dg(z'), sup/ (d M {z',z")y 2 dg(z") < Ce 

z" J d M (z' ,z")<e z' J d M (z' ,z")<s 

so by Schur's test we have 

\u(t,z')\ 2 \u(t,z")\ 2 (d M (z',z"))- 2 dg(z')dg(z") 

(z',z")<£ 

<Ce\\\u(t)\ 2 \\ LHM) \\\u(t)\ 2 \\ LHM) . 

Integrating this in t we see that this contribution to l|4.9|l is 0(e J M \u(t, z)\ 4 dgdt) 
= O(Z) as desired. 

Now we prove (|4.10() . From i|3.3[) we have 



\u{t,z")\ 2 dg{z") = O{\\u{0)\\h { M)) 
so we can bound the left-hand side of (|4 . 1 01) as 



O(\\u(0)\\l HM) / sup I / 1 '^ ( \: 2 ' 3 Vl ' dg(z')\dt) 
Jo z"eM Jm [ r ) 

which by Lemma 14.31 (letting the index set W be the manifold M) and l|3.3l) is 

O(|| U (0)||| 2(M) || U (0)||^ 1/2(M) ) = O(2 k \\u(0)\\t 2(M) ) = O(Z) as desired. 

Finally we prove (|4.11() . From (|3.3|l . (|3.4() and Cauchy-Schwarz we have 

\u(t,z")\\Vu(t,z")\ dg(z") = O(|| U (0)|| i2(M) || M (0)|| H i (M) ) = O(2 fe || U (0)|| 2 L2(M) ) 

so by using Lemma |4~31 and (|3.3(l as with (|4.10() we can bound the left-hand side of 
gUll as 

o(2 k \\umi H M)\wmh { M)) = o(z) 

as desired. □ 

5. i 4 ESTIMATE IN THE NEAR REGION 

In this section we prove the near region estimate Proposition 14. II which is sig- 
nificantly easier to prove (and indeed also follows from earlier work in [201, [2]) but 
already illustrates the basic method. We fix rj and zo, and allow our constants to 
depend implicitly on zq, rj. 

Let A n ear be the commutant defined in (|4.3(l , and let U — u <E> u as above. From 
(12.61) . Cauchy-Schwarz, (|3.3|) . (|3.4|) and the trivial observation that c?m is Lipschitz, 
we see that 

I (Aicar U{t),U(t)) MxM \ < C\\U{t)\\ H l( MxM )\\U(t)\\ L 2( M> , M ) 

<C2 k \\u{Q)\\U M )=0(Z) 
for t = 0, 1. Thus by (|2~T|) we have 

(5.1) / (i[H M xM, A aeai ]U{t), U(t)) M *M dt < O(Z). 



The idea is now to decompose the left-hand side of 15.1|l as a positive main term 
that will yield Proposition 14. II plus some error terms which are either positive or 
0(Z). 
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We first expand the commutator using (|4.3|l as 

near J — [Hm xm, <p][Hmxm, d>M( z ', z ")] 
-<p[H M xM, [Hm xm, Am{z\ z")]]. 

Consider the contribution of the first term [HmxM, ¥>][Hmxm, d,M(z', z")] in (|5.2() 
to ()5.1(l . This is an error term and will be bounded in absolute value. By integration 
by parts, we can bound this contribution by 

(5.3) 0([ \([H Mx M,d M }U{t) 7 [H Mx M^}U{t))MxM\dt). 

Jo 

To control this expression we apply 1(2.2(1 to expand 
(5.4) 

[H M xM,d M (z',z")] = (V z >d M ,Vz>)g + ((V z »d M ),V 2 ») s + (H z ,d M ) + (H z „d M ). 

To estimate these quantities, we first recall the form of the Laplacian in local polar 
co-ordinates: 

Lemma 5.1. Let x % be normal co-ordinates on a Riemannian manifold X, cen- 
tred at p € X, and let (s,(f>) be polar co-ordinates with respect to x l . That is, 
s = X^C 2 - 4 ) 2 an d are homogeneous of degree zero with respect to x 1 . Then the 
Laplacian takes the form 

d 2 .n — ld 2 . 

os z s os 

Here c is a smooth function which is 0(s 2 ) and A s »-i is the Laplacian on S"^ 1 

determined by the co-ordinates (j> and the metric s^ 2 g restricted to the submanifold 
{s = s }. 

Proof. The Laplacian in Riemannian polar normal co-ordinates is given by 

— dsyfgds + -^pd^iyfgg^d^, g = ds 2 + h 3 (s, <j))d(j) 1 dep 3 , h 3 = 0(s 2 ). 
y9 v 9 

Since in normal co-ordinates the metric is Euclidean to second order at the origin, 
^fg = s™ _1 (l + 0(s 2 )). The result follows readily from these facts. □ 

From this and 15.4J1 we see that in the region z', z" G B(zq, -q) we can write 

i[H M xM,d M (z',z")} = 0{l)V z ,, z „+0{l)/d M {z',z") 

where 0(1 ) denotes various bounded functions of z', z" . Applying this estimate, as 
well as using (12.2(1 to expand i[HmxMi <p], we can thus bound ((5.3(1 by 

0{C [ (\U(t,z',z")\ + \W z/ ^U(t,z',z")\) 

JO J B(z a ,2ri)xB(z„,2ri) 

(\U(t, z', z")\/d M (z', z") + \V z ,, z „U(t, z', z")\) dg(z')dg(z")dt.) 

Expanding out U and using Cauchy-Schwarz, we can bound this by a linear com- 
bination of (|4.7(l . 14.9(1 . Thus these terms are O(Z). 

Comparing the above estimates with 1(5. l|l and 15.2(1 . we now have 

(5.5) f (-<p[H MX M, [H M xM,d M }}U(t), U{t)) M xM dt < O(Z). 

Jo 
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Next, we apply l|2.3(l to obtain 

(5.6) -[H MxM , [H M xM,d M }} = -V Uess(d M ) af) V a ~ H 2 MxM d M 

where the gradients V are on the product manifold M X M (so the indices a, (3 
range over six values). 

Now consider the contribution of the lower order term H 2 , IxM dM{z' , z" ) to l|B~5|l . 
We need 

Lemma 5.2. For z' \z" £ B(zq,2i]), we have 

~H 2 MxM d M (z',z") = 8ir6 z ,(z") + O^ 1 ). 

Proof. We again use polar normal co-ordinates. Choose a co-ordinate system so 
that for each fixed z" the z 1 co-ordinates are normal centred at z" . Then changing 
to polar co-ordinates s, <j> in the z' variable we see that 

A z ,d M (z>, z") = (-* - (l + c)^4- + s- 2 A s ,-A s 



ds 2 s ds 

, n — 1 
= 1 + c , 

s 

where c(z! , z") vanishes to second order at z' = z". This condition is symmetric 
under interchanging z' and z" so the same is true of Amx-M^m^', z"); specializing 
to n = 3, we thus have 

A MxM d M (z', z") = (4 + c)d M {z', z")- 1 

where c also vanishes to second order at z' — z" . Now we apply Amxm again using 
Lemma f5. II to obtain 

AMxMA M xMdM(z',z") = -32n5 z ,(z") + 0(d M (z' , z")- 1 ) 

and the claim follows (note that this is consistent with l|2.9[l ). □ 

Thus the contribution of this term to (|5.5|) is of the form 

8ir f [ <f{z,z)\u(t,z)\ 4 dg{z) 
Jo Jm 
"i r 



+ 0( / \u(z")\ 2 \u(z')\ 2 /dM(z,z')dg(z')dg(z")). 

JQ J B(z ,2-q)xB(z (u 2ri) 

The first term is the main term (as in (|2.9|) and ^.lOp l. The error term is 0{Z) by 
(Oil . 

The remaining term to consider in l|5.5|l is the contribution of — VHess(<i/vf )V, 
which we can write using ljl.5|l and integration by parts as 



(5.7) 



Re / / pE.ess(d M ) af} V a U(t)VpU(t) dg(z')dg(z")dt 

JO JMxM 



Re / (V / g^)Hess(d M ) Q/3 V Q ?7(i);7(i) dg(z')dg(z")dt. 

JO JmxM 



The first term in 1)5. 7|l is the most interesting one. We use the characterization 
of Ress(dM)(v, v), where v = (v',v") £ Ti z i , Z »)M 2 , as the second derivative of 
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dM{z'(t), z"(t)) where (z'(t), z"(t)) moves along a geodesic in M 2 with initial con- 
dition ((z', z")\ v). We use the second variation formula for geodesies (see |12j . 
Theorem 4.1.1) to obtain 

(5. 8) P«w<m > -35^5 jf {m,#v*.<>.)« 

> -c<i M ( a '(t),/'(t))(i»T + i«"i a ), 

where d s is the velocity of the geodesic which goes from z'(0) to z"(0) in unit time, 
hence |<9 S | = dn/j(z', z") and dt is the Jacobi field corresponding to the variation of 
the geodesic. In particular this bound does not blow up when di^iz 1 ,z") is small. 
From this it follows that 

Re / [ ^Hess(rfM) Q/3 V a C/(t)V U{t) dg(z')dg(z")dt 

J0 J M x M 

> -C [ f |VC/(t)| 2 dg(z')dg(z")dt. 

JO J B(za,2r])xB(z ,2r]) 

The right-hand side is then O(Z) by g2J. 

Finally consider the second term in l|5.7|l . Using the crude estimate |Hess(dA/)| < 
C I du we can bound this contribution by 

0{ f f \U(t,z',z")\\V z ,, z „U(t,z',z")\/d M (z',z") dg(z')dg(z")dt). 

Jo J B{z„,2ri)xB(zo,2ri) 

By Cauchy-Schwarz we can bound this a linear combination of (|4.7|l and (|4.9|l , and 
so these terms are also O(Z). 

Inserting all of the above estimates into Q5.5|) . we obtain Proposition ^. II 

6. The geometry of the scattering region 

To prove Theorem 11.51 it remains to verify the scattering region L 4 estimate 
in Proposition 14.21 To do this we shall need to understand the geometry of the 
scattering region. More precisely, we will need to control the metric function 
du{z' ', z") — du((x' ,y'), (x",y")) and its derivatives on the support of \ip, where 
X, ip are as i n Q4.4[l . To avoid singularities forming in dM, we will always assume 
in this section that tq is sufficiently large and rj sufficiently small. 

We begin with a basic symbol estimate on the metric dM- 

Proposition 6.1. Let (z',z") be in the support ofxi>- If du[z' , z") > - +q , then 
we have the product symbol estimates 

\^ivi,\ z >, z nd M {z',z")\ < C m ,, m/ ,(r')- m '(r")- m " 

for all m! \m" > 0. If instead duiz' , z") < r ^ , then we have 

\Vy^,d M (z',z")\ < C m d u {z',z") 1 - m 

for all m > 0. Furthermore, for j < 3, the distribution V^, z ,,dM{z / 1 z") has no 
mass on the diagonal z' = z" . 

This symbol estimate, while quite plausible, is a little technical and its proof will 
be deferred to Section El From this, Lemma 15.21 and a rescaling argument we can 
obtain the following estimate for the bi-Laplacian of dM- 



20 



ANDREW HASSELL, TERENCE TAO, AND JARED WUNSCH 



Corollary 6.2. The function Hj^ xM d,M{z' , z") satisfies the distributional inequal- 
ity 

-H 2 IxM d M (z',z") < 8n6 z ,(z") - Cd M (z', z")"V + r")" 2 - C(rT 3 - C(r'T 3 
uniformly on the support ofx4>- 

Proof. Observe that if \r"/r' — 1| > c for some c > (Im{z',z") is comparable 
to r' + r" and the estimate follows from Proposition 16. II so we may assume that 
\r"/r' — 1| < c for some small c. In this region Proposition 16. II is insufficient (it 
gives a bound of 0(djf), which is not even locally integrable near the diagonal). 
Instead, we we set t := l/r' and consider the rescaled manifold 12 M T with metric 

g T = dr 2 + r 2 h ij (- 7 y)dy t dy J . 
r 

Observe that g T various smoothly in r as r varies in the compact set [0, 2/ro]. Also, 
since ddM(y',y") < 2?? and \r"/r' — 1| < c, we see that 

d MT {Tz\rz") < C(c + n), 

with the action tz' := (rr',y') and likewise in z" . Thus if c and r\ are sufficiently 
small, we can apply Lemma 15.21 and conclude that 

-H 2 MTXMT dM r (z',z") < M5 tz ,{tz") - CduArz'^rz'T 1 

uniformly in r. The claim follows by undoing the scaling. □ 

Finally, we need a geodesic convexity estimate on the metric dM(z'(s), z"(s)), 
generalizing Q2.8JI . 

Proposition 6.3. Let z'(s), z"(s) be two geodesies. For ro sufficiently large and rj 
sufficiently small, the estimate 

d 2 Av' I 2 |v" l 2 \ 

^d M (z'(s), z"{s)) > -C(i^- + !^-) 

(6 ' 1} r /KI 2 , \v'T\ 

V(r') 2 {r") 2 J 

holds whenever z'(s),z"(s) are distinct points in the support of xty, where \v'\ := 
\£ z '( s )\g(z'(s)) is the speed of z' and |w^ ng | := r'(s)\-£y'(s)\ h ^^ is the angular 
component of this speed, and similarly for v" . 

This proposition is somewhat delicate and will be also deferred to Sectional The 
idea is to first prove this lemma for perfectly conic metrics (in which the second 
group of terms in l|6.1|l do not appear) ; the perfectly conic metric contains a large 
number of totally geodesic planes on which the metric behaves like the Euclidean 
metric except in normal directions (generating the first group of terms in 
We then handle the approximately conic case by obtaining symbol estimates for 
the error between the approximately conic and perfectly conic metrics giving rise 
to the second group of terms in (|6.1|l . 



Strictly speaking, this rescaling is only well-defined on the scattering region of the manifold, 
but that is all we shall need here. 
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7. L 4 ESTIMATE IN THE SCATTERING REGION 

We can now begin the proof of Proposition 14.21 Fix ro,rj; we assume ro to 
be large enough and 77 small enough that the geometrical lemmas of the previous 
section apply, and allow our constants C to depend on ro, 77. Let A asymp t be as in 
1)4.4(1 . By arguing as in Section 03 we have 

(7.1) f (i[H M xM,A asympt }U{t),U(t)) MxM dt<0(Z). 
Jo 

Again, we will expand the left-hand side of 1(7.1(1 to extract a positive main term 
that will give the estimate in Proposition ^. 21 plus a number of error terms which 
are either positive or O(Z). 

From l(4.4|l we may decompose 

i [Hm x M , ^asy mpt ] = — [H M x M , X] i> [H M x M j ] 

(7.2) -x[HMxM,ip][HMxM,d M ] 

— Xl/)[i?MxM, [ffMxM,4]]- 

The first two terms are error terms and will be estimated in absolute value; it is 
the third one which will generate the interesting positive terms. 

We first analyze the first term [Hmxm , xI^I^Mxm , du] in 1)7.2(1 . Integrating by 
parts we see that the contribution of this term to 1(7.1(1 is 

(7.3) 0([ \mHMxM,x}U(t),i[H MxM ,d M }U(t)) MxM \dt). 

Jo 

This integral is supported in the region where either (z 1 ) < ro or (z") < ro, since 
\ is constant elsewhere; without loss of generality it suffices to consider the terms 
when (z') < tq. To control this expression, we first observe from 12.2(1 and Lemma 
IO that 

(7.4) i[HMxM,d M ]=0{\)V z ,, z n+0{\)/d M (z',z") + 0{l) 

on the support of x^i where we use O(l) to denote various bounded (tensor-valued) 
functions of z', z"; for future reference we also remark that these O(l) errors obey 
symbol estimates when ddM (y', v") > »?/2, since in this case d(z', z") > c v (r' + r"). 

From ((7.4(1 . 1(2.2(1 and Cauchy-Schwarz, we thus see that ((7.3(1 can be bounded 
in magnitude by 

0( [ f \Vu(t, z')\ 2 \u(t, z")\ 2 + Ht, z')| 2 |V«(t, z")\ 2 

Jo J{z')<r 

+ |u(i, z')\ 2 \u(t, z")| 2 (l + d M (z', z")- 2 ) dg(z')dg(z")). 

The first two terms are O(Z) by ((4.7(1 . For the third term, observe that the con- 
tribution of the region r" < 4ro is O(Z) by 14.9(1 . while for the region r" > 4tq we 
simply estimate (1 + dm(z' , z")~ 2 ) crudely by O(l), and this term is then bounded 
using ((3.3(1 by O(||m(0)||^ 2 (a/)) = O(Z). Thus the total contribution of the first 
term of (72) to £7JTJ is O(Z). 

Now consider the contribution of the second summand of 1(7.2(1 to ((7.1|) : these 
are also error terms, but of a more delicate nature. Again we integrate by parts to 
estimate this contribution by 

(7.5) 0(f \(i[H Mx M^}xU(t),i[HMxM,d M }U(t))MxM\dt). 

Jo 
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The key point here is that since the angular cutoff ip does not depend on the 
radial variables r',r", the first-order operator [Hmxm,iP] only contains angular 
derivatives V y ^ y " and no radial derivatives V r ' )r », thanks to l|2.2|l . Indeed, on the 
support of x we have 

[H m ,m^] =0(l)^V^ + 0(l)ppV^ 

where again we use 0(1) to denote various bounded functions of z' , z" which also 
obey symbol estimates. From this, <|7.4[1 . and symmetry we see that these contri- 
butions are bounded by a combination of expressions of the form 



, , , a^ k (z',z")fXt^')V k u(t,z') „ , , „ ,. 

oil — ^ u t)Z //)|2 dg{z >) dg{z ») dt , 



JMxM 



a? k {z' ) z")f i u(t,z r )u(t,z') 



0(\ / ^ T — -u(t,z")V k u{t,z")dg(z')dg(z") dt\ ), 

V Jo Jmxm r 

of | f I *|Y 

V JO Jr'.r">r r'd M {z',z") 'J 

(\ f [ |V»(M0H^)ll»(M")l 2 , t| \ 

V JO Jr',r">r v) ) 
f, f 1 f \u{t,z')\ 2 \u{t,z")\ 2 ,M ,| 

(noting that V a = rj7), where a^ k denotes symbols supported on the the support 
of x^V' (note that we do not encounter the singularity on the diagonal z' = z" 
because V?/> vanishes here). The first two terms are controlled by (|4.1U|I . (I4.11|) 
respectively. The third and fourth terms can be controlled via Cauchy-Schwarz by 
a linear combination of (|4.7|l and (|4.9|) . Finally, the last term is easily controlled 
by (|4.9|) . (Recall that r' is large.) Thus the total contribution of the second term 
of JE21 to (□} is O(Z). 

In light of H7.1(l . Ij7.2|) . and the preceding estimates, we see that 

(7.6) f (- X iP[Hmxm, [H m *m, d M (z', z")]]U{t), U(t)) MxM dt < O(Z). 
Jo 

Of course, we may expand — [Hmxm, [Hmxm, dM(z' , z")}] using l|2.3|l as 

(7.7) -{H MxM , [H MxM ,d M (z', z")]] = -V^Hess(d M ) Q/3 V Q - H 2 MxM d M 

where Hess(d M ) Q/3 = VV^d. 

Let us first study the lower order term —H\ IxM dM{z',z"). Applying Corollary 
16.21 the contribution of this term to 17.611 is greater than or equal to 

> 8tt / / X (z,z) 2 \u(t,z)\ i dg{z)dt 
/o Jm 

\u(t,z')\ 2 \u(t,z")\ 2 , , ,, , , „,_, , 
1 V, , / ' L dg(z' dg(z")dt). 

The error term here is of course O(Z) by (|4.9|) . 



+o< /7 

JO JM 
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It remains to treat the contribution of the — VHess(dM)V to Ij7.6|) . Again we 
integrate by parts, creating a main term 

(7.8) Re f f (x(z, z")i/>(y', y")Uess{d M ) a0 V a U{t, z' , z")V p U(t, z', z") dt 

JO J MxM 

and an error term which is bounded by 

/ |V z%z »(x(^^'')V'(2/',/)l|Hess(dAf)||V^^!7(i,z',z'')||C/(t ; z',z'')l dt. 

MxM 

Consider the error term first. For this term we use the crude estimate 

\Ress(d M )(z' , z")\ < C(l + d M (z, z")- 1 ) 
from Lemma l6.ll as well as 

\V z ,, z „(x(z l ,z")^(y , ,y")\<^ + ^- 

By symmetry we can thus estimate this error term by 

C f [ \Vu(t,z')\\u(t,z')\\u(t,z'')\ 2 (^ + ^)(l+ \ ) dg{z')dg{z")dt. 

JO J MxM r r Um(Z,Z ) 

By Cauchy-Schwarz one can estimate this by a linear combination of (|4.7|l and 
(Oil . 

Now consider the main term (|7.8I) . Recall that Hess(djvf)(u, v), where v <G 
T( z / jZ //)M 2 , is the second derivative of (Im(z'(s), z"(s)) in s as (z'(s), z"(s)) moves 
along a geodesic in M 2 with initial condition ((z' , z");v). By Proposition 16.31 we 
see that 

f\fu(t,z>)\i\u(t,z")\ 2 , \u(t,z>)\ 2 \fu(t,z")\ 2 \ 



(EHJ > - c 



r ■ ,r >ro /2 

V«(t, z')\ 2 \u(t, z")\ 2 \u(t, z')\ 2 \\7u{t 7 z")\ 2 



dg(z')dg(z")dt. 



\ (r 1 ) 2 ^ r ")2 

But the first term on the right-hand side is O(Z) by l|4.8|l . while the second term is 
also 0{Z) by l|4.7() . Combining all of these estimates together we obtain Proposition 
14.21 (once we verify Proposition 19.41 and Proposition 16.31 in the next section) . This 
completes the proof of Theorem 11.51 

8. Concluding remarks 

We first want to comment on some particular geometric settings in which the 
proof of the Morawetz estimate is much easier. The first is that of a simply- 
connected smooth three-manifold M which has globally non-positive, bounded sec- 
tional curvature; by the Hadamard-Cartan theorem (see ^21> Corollary 4.8.1), such 
a manifold is diffeomorphic to R 3 . Such a manifold is automatically non-trapping, 
and the distance function du{z\ z") is smooth on all of M x M outside the diago- 
nal. We also make the technical assumption that the distance function du{z' ', z") 
is uniformly C 4 in the off-diagonal region {(z',z") G M x M : du{z' ' , z") > l} 13 . 
In particular H 2 dM is bounded in this region. More importantly the distance func- 
tion is globally convex (see ^2]; Theorem 4.1.1), so that the tensor Hess(dM) al3 is 
positive semidefinite, by the second variation formula. Thus with the commutant 



13 It may be that this is automatically true, but this is not clear to us. 
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A = i[H, cIm(z') z")\ (with no cutoffs whatsoever) we obtain a Morawetz inequality 
very easily for this class of manifolds. 

The next case we consider is that of asymptotically flat three-dimensional non- 
trapping manifolds, i.e. where M = R 3 and the metric g satisfies with h equal 
to the standard metric on S 2 . We make a distinction between the distance d(z', z") 
on M induced by the metric g and the Euclidean distance \z' — z"\ obtained by 
identifying M with Euclidean space. We shall only discuss the Morawetz estimate 
in the scattering region (z'), (z") > 7*0. Here the key observation is that when z 1 
and z" are fairly close together, e.g. \z' — z"\ < ^\z'\, then the distance duiz! ,z") 
is smooth and is close to \z' — z"\ in the sense that 

\cIm (z',z") — \z — z"\\ < C whenever (z 1 ), (z") > tq and \z' — z"\ < ^l 2 '!' 

at least if we make the compact region Kq sufficiently large. This is an easy con- 
sequence of the decay of the difference g — 8 between the two metrics. In fact, the 
error obeys symbol estimates of one order better than Proposition lfi. II suggests: see 
Proposition 19.41 in the appendix. We now consider the commutant A := i[H,a], 
where a is the function 

(8.1) a := x(z', z"){ip{z' , z")d M (z', z") + (1 - ip{z' , z"))\z' - z"\), 

if is a smooth cutoff which equals 1 when \z' — z"\ < j\z'\ and equals when 
\z' — z"\ > h\z'\, and x i s as m Sectional Thus this cutoff is equal to the actual 
distance d(z', z") when z' and z" are close, but reverts smoothly to the Euclidean 
distance when z' and z" are far apart. This significantly reduces the error caused 
by differentiating the cutoff cp, indeed it generates terms that are now controlled 
by l|4.7(l rather than the terms l|4.9jl . 1)4. !()[! . (|4. 1 1|) which required the Morawetz 
estimate ()l().6|l and its refinement in Lemma f4. HI We omit the details. 

Next, we give an example to show that 14.5fl cannot be strengthened by bringing 
the absolute value signs inside the ^-integral. 

Proposition 8.1. There does not exist an estimate of the form 

, ^ f 1 f \fu(t,z)\ \Vu(t,z)\ „ nh9 

(8.2) / / Lg v ^11 \±Jl dgdt < C \\u(Q)f H1/Hm (FALSE!) 

JO JM r 

for solutions to (|1.2[> even in Euclidean space (M,g) — (R 3 ,5). 

Proof of falsity of 1)8.2(1 . We give a sketch only. Pick a large integer N = 2 k > 1, 
and let <j> be the explicit solution to 1)1.2(1 given by 

4>(t,z) := (t + i)- 3/2 e^^; 

note that at time t = this is basically a standard Gaussian, while for times t ^> 1, 
the solution <j>(t) is concentrated on a ball of radius 0(t) centered at the origin, and 
it (and its derivatives) have magnitude roughly t~ 3/>2 on this ball. We then let u 
be the function 

fe 

u(t, z) := k^^NH, Nz) + e j( j){NH, N(z - 2 j ei )) 

i=i 

where e\ = (1, 0, 0) is a standard unit vector and ej = ±1 are random signs chosen 
independently; the purpose of the signs is to ensure (thanks to Khinchin's inequal- 
ity) that we do not expect any unusual cancellation between the summands. 
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The function u{t,z) is clearly a solution to H1.2(l . and the ff 1 / 2 (R 3 ) norm can 
be computed as 

IKOJII^/^-fcVajvr 1 ; 
this is easiest shown by first verifying the more general statement || , u(0)||jj»(r.3) ~ 
A; 1 / 2 .^ -3 / 2 for s = 0, 1, 2 (applying a rescaling by N if desired) and then interpo- 
lating. Thus the right-hand side of (JOl is 0{kN- 2 ). 

Let us choose any iV/3 < j < 2N/3 and consider the size of u(t,z) and its 
derivatives in the region of spacetime where t ~ 2 J /A r and |z — 2 J ei | <C 2 J . In 
this region |V/c 1 / 2 0(iV 2 i, iVz)| has size ~ fc 1 / 2 (2 3 N)~ 3 / 2 , and so the expectation of 
|Vu(t, z)| is also at least > k 1 ' 2 (2 J N)~ 3 ' 2 . Since <fi is radial, the angular derivative 
jyfc 1 / 2 ( /)(7V 2 t, Nz)\ vanishes, however \f(f){N 2 t, N(z - 2 1 e 1 ))\ is fairly large, com- 
parable to ~ (2 : 'N)~ 3 ' 2 in a large fraction of the region under consideration. Thus 
this region of spacetime contributes at least > k 1/2 (2 j N)- 3 (2 j N)2 3j /2 J ~ fc^jV" 2 
to the integral in l|8.2J) . Summing over all j we see that the left-hand side of (|8.2J) is 
at least ~ fc 3//2 7V -2 . Comparing this with the right-hand side we obtain the desired 
contradiction 1 by setting k — > oo. □ 

Remark 8.2. One can construct a similar example to show that the second estimate 
in Lemma 14.31 similarly fails if the absolute values are placed inside the integral. 

Remark 8.3. Note that in our example y7u(t,x) and (the largest term in) \7u(t,x) 
oscillate in different directions and so if one removes the absolute values, replacing 
instead by a smooth symbol, then we do not obtain a counterexample to Lemma 

m 

Finally, we remark on the applicability of the Strichartz estimate in Theorem ll.5l 
to non-linear Schrodinger equations. It turns out that in such equations, the most 
useful Strichartz estimates are either those that require no derivatives whatsoever 
on u (i.e. they have L 2 (M) on the right-hand side), or measure a in a space of the 
form L q t U^ . Our current estimate has neither of these two properties. However 
by commuting the Schrodinger flow with a power of (1 + H) and using Sobolev 
embedding one can get an LfL^ estimate, for instance we have 

IMIl4l~([0,1]xM) < C £ ||u(0)|| ff i+ E ( M ) 

for any e > 0. This result can for instance be used to demonstrate local well- 
posedness in H 1+e (M) of the quintic non-linear Schrodinger equation iut + ^Au = 
±|u| 4 u on non-trapping asymptotically conic manifolds by a standard argument 
which we omit (see e.g. jS], |3])- Note that the above estimate was also derived by 
Burq |2] for asymptotically flat non-trapping manifolds. 

It would be interesting to see if the e can be removed; this could then be used 
to demonstrate global well-posedness of the above quintic equation in the energy 
space H 1 (M) for small energy data, using the conservation of energy in the standard 
manner (note the sign of the non-linearity is irrelevant in the small energy setting) . 
While the e can indeed be removed from the above Strichartz estimate in the 
Euclidean space case (see [21]; the corollary concerning global well-posedness can 
be achieved by many other means, see for instance HJ), we were unable to remove 

14 Note that this is only a 'logarithmic' failure in the frequency variable. Indeed if we concede 
an epsilon worth of derivatives then we can handle these terms easily by a variant of 11.111 where 
the epsilon loss is transferred from the weight r~ 1— e to the regularity i? 1 ^ 2 + e (M), thus recovering 
the result of Burq 0, at least in the context of Lf x estimates. 



26 



ANDREW HASSELL, TERENCE TAO, AND .TARED WUNSCH 



it here. However, by a Besov space interpolation argument (see [23) it is possible 
to obtain the variant estimate 

ll M IU?L~([0,l]xM) < C«||w(0)||h 3 / 2 - 2 /(9-i)(M) 

for all 4 < q < oo. This suffices to give local well-posedness on non-linear Schrodinger 
equations iu t + \ Au = F(u) in the space H 3 / 2 ~ 2 / (M) for any q > 4, whenever 
F(u) is smooth and grows like \u\ q ; we omit the details as they are rather standard 
(see e.g. [5])- Note that in Euclidean space the space 

#3/2-2/(9-1) WQuld be criti _ 

cal with respect to the usual scaling of the non-linear Schrodinger equation. Thus 
the Strichartz estimate Theorem ll.5l has some application to non-linear Schrodinger 
equations, though clearly it is less satisfactory than the full range ljl.8|l of Strichartz 
estimates would be; these would yield local well-posedness in H s for any q > 1 and 
s > max(0, § - ^-). 



9. Appendix I: Proof of geometrical lemmas 

In this section we analyze the behaviour of geodesies in thin conic regions of 
M, in order to prove the geometric estimates in Proposition 16.11 and Proposition 
16.31 required in the above proof of Theorem 11.51 The main strategy is to compare 
the asymptotically conic metric cLm to a perfectly conic metric d con ic, which has an 
exact formula and can be handled explicitly. 

By compactness of dM it suffices to work in a truncated sector 

(9.1) &(yo,ri> r o) : = {( r , y) I d dM {y, yo) < m ( z ) > M 

for some yo S dM, where we assume ro sufficiently large and r\ sufficiently small 
(possibly depending on yo, although this is irrelevant since yo ranges in a compact 
set). 

We first briefly review geodesic flow on a general Riemannian manifold (M, g) 
parameterized by a general set of co-ordinates zK We will view geodesic flow as 
a Hamiltonian flow on the cotangent bundle T*(M) of M, which we parameterize 
as (z^Xj), where Q is the co-ordinate dual to the vector field dz^; thus Qdz? is 
the canonical one-form. On this cotangent bundle we define the energy function 
a(H) : T*(M) — > R as 

a{H) :=^ fe (z)Oa = i|CI^); 

note that this is also the symbol of the Schrodinger operator H. This energy 
function then induces trajectories (z(s), £(•?)) in the cotangent bundle which are 
just the paths of geodesic flow (using the metric g to identify the cotangent and 
tangent bundles): 

^ - -^ex*- 

These flows preserve cr(H), and are also homogeneous with respect to the scaling 
(z(s), C( s )) l— * {z{^s), X((Xs)), which corresponds to speeding up the velocity of the 
geodesic flow by A. If we use this homogeneity to normalize to the unit speed case 
a(H) = i (i.e. ICIgO) = 1)) then this flow is related to the distance function dM on 
M by the formula 

d M {z{s),z(s')) = \s- s'\ 
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provided that \s — s'\ is sufficiently small. Thus the short-time geodesic flow can 
be used to define the metric locally. 

Now we specialize to studying geodesic flow in the truncated sector £l(yo, n, ro) in 
the scattering co-ordinates (x,y J ), where the manifold M is either asymptotically 
conic or perfectly conic. Of course this region is not geodesically complete, so the 
statements below are conditional assuming that the geodesic does indeed stay inside 
this region. 

The first step is to choose co-ordinates for the cotangent bundle T*f2(?/o, Vi r o)- 
We shall use the co-ordinates (x, yi , v, \i.f), where < x < e , and y g dM, while 
the co-ordinates v and [ij are dual to the vector fields —x 2 d x — d r and xd y j 
respectively 15 ; thus one can think of n as an element of T*(dM), the canonical 
one-form is —vdx/x 2 + ^dyj/x. Following Melrose we define the scattering 
cotangent bundle over the compact manifold M as the bundle SC T*M whose sections 
are locally spanned over C°°(M) by dx/x 2 and dy/x (hence can be paired with the 
scattering vector fields, spanned by —x 2 d x and xd y ). Because of the form of the 
metric, the symbol cr(H) of the Schrodinger operator H now has the form 

°{h) = \{v 2 + m*)mV) = + H£ W ) 

and thus in the unit speed case o~(H) = 5, v and \i will have combined magnitude 
equal to 1. 

The geodesic flow on SC T*M can then be written explicitly as 
d 2 

— X = —x V 

ds 

d A . jU 

— y J = xh J \x k 

(9 ' 2 ^ d L , fc 1 2 dh? k 

— v = xh J fijfik + 7^ "faTVjVk 

d 1 dh kl 

recall that h? (x, y) is, for small x, a family of metrics in y with x a smooth 
parameter; it is only in the perfectly conic case that there is no ir-dependence. 

This geodesic flow of course still preserves <j(H), and is still homogeneous with 
respect to the scaling (x(s), y(s), v(s), //(s)) 1— ► (x(Xs), y(As), Xv(Xs), A/i(As)). This 
geodesic flow equation can be used to define the distance dudx', y'), (x", y")) pro- 
vided (x',y r ) and (x",y") are sufficiently close; we will quantify what 'sufficiently 
close' means later. 

To gain some intuition as to how the distance du behaves on Q(yo, n, ro), let us 
first consider the case hP k (x,y) = h^ k {y) when M is perfectly conic near infinity; 
to distinguish this from the asymptotically conic situation we shall write (iconic = 
iconic, h instead of c?m, and write the metric g as iconic. h- In this case the equations 
for geodesic flow are equivalent to those on a plane, for the following geometric 
reason. Consider a geodesic 7 in BQ^j(y Q ,ri), and consider the truncated cone 
{(x,y) : < x < e ,y E 7} over that geodesic. In the perfectly conic case, this 



15 More informally, a symbol such as a(x,y,u,/i) corresponds in the usual polar co-ordinates 
r, y to the pseudo-differential operator a(i, y, J^;,)P). The advantage of these co-ordinates is that 
v, p, remain bounded for unit speed geodesic flows. 
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truncated cone is both perfectly flat and totally geodesic , and in particular is 
isometric to a subset of the plane. Conversely, it is easy to show that there is a 
unique minimal geodesic connecting (a/, y') and {x" , y") which lies in the truncated 
cone over the geodesic in dM connecting y' and y" . From plane geometry we thus 
see that we have the cosine rule 17 

(9.3) d conic (z', z"f = (r') 2 + (r") 2 - 2rV" cosd aM (y', y") 

for z',z" € Q(yo,r],ro). Note in particular that we thus expect symbol-type regu- 
larity estimates on clonic, in the sense that each derivative in the z' and z" variable 
gains a power of d; this is of course what happens in the Euclidean case, and we 
will be able to show that it also happens in asymptotically conic manifolds near 
infinity. This is quite plausible when r' and r" are of comparable magnitude; the 
most delicate issue will be proving these sorts of bounds when r' and r" are both 
large but far apart from each other, which seems to be a geometric necessity in 
order to not lose any derivatives in Theorem 1 1.51 

This completes our discussion of the perfectly conic case. We now consider the 
distance function dM of an asymptotically conic manifold on fi(j/oj Vi r o)- We shall 
always take rj to be sufficiently small and r$ to be sufficiently large. 

We first show that globally geodesies between points in Q(yo, r], ro) stay in the 
scattering region. 

Lemma 9.1. If rj is sufficiently small and tq sufficiently large, then for every pair 
of points z' ,z" in fl(jjo,rj,ro), every globally minimizing geodesic between them lies 
entirely inside £l(yo, 5ry, ro/2). 

Proof. Without loss of generality we may take r" > r' . Suppose for contradiction 
that the globally minimizing curve connecting z' and z" contains a point (r, y) in 
the complement of Cl(5r], r o/2). Then either r < ro/2 or dgM(y,yo) > 5rj- 

First suppose that r < ro/2. Using the crude estimate g > dr 2 , we thus obtain 
the lower bound 

(9.4) d M (z\ z") > (/ - r /2) + (r" - r„/2) > r". 

On the other hand, if ro is large enough then we have the pointwise bound hij(l/ r, y) 
< ihij(0, y) for all r > r /10, y G 9M, and hence g < g CO nic,4h m ths region. Since 
the minimal geodesic connecting z' to z" in the iconic, Ah metric lies in the region 
r > r o/10, we thus obtain from 1)9.3(1 

(9.5) d M (z', z") < d conicAh (z', z") < v /(r') 2 + (r") 2 -2rV'cos4 ?? , 

which contradicts (|9.4|l if n is sufficiently small. 

Next suppose that ddM{y,yo) > 5??. Then the distance from y to both y' and 
y" must be at least 4ry. Using the comparison /iy(0, y)/4 < hij(l/r,y) < 4hij(Q,y), 
which holds for r > ro/2 if ro is sufficiently large, we have the pointwise bound 



Indeed, one can think of 7 as being isometric to an arc in S 1 , so the truncated cone is then 
isometric to a truncated sector in R. 2 . 

^This rule can also be deduced from 19.21 : the point is that the term in the equation 

for -^v disappears in the perfectly conic case, and so the geodesic flow equation decouples into 
a (rescaled) geodesic flow equation on dM for y,fi, and an explicitly solvable evolution for x, v. 
We omit the details. 
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dr 2 + hij (0, y)/4 < g < dr 2 + Ahij (0,y). Comparison with this smaller conic metric 
shows that the length of the curve must be at least 

(9.6) V( r ') 2 +r 2 - 2r'r cos 2n + v / {r 2 + (r") 2 - 2rr" cos 2rj, 

since the broken geodesic in the small conic metric that goes between these points 
has at least this length. For r) < tt/12 this is larger than 1|9.5|) . as can be seen by 
consideration of the lengths of a plane Euclidean triangle with vertices having polar 
coordinates (r',0), (r",4ry) and (r',2rj). This is a contradiction. □ 

We now use a symplectic argument to deduce regularity of the distance function 

cIm(z', z") where (z 1 , z") lie in some cone Q(y, rj, ro) for fixed (77, ro). To do this, we 

fully exploit the smoothness of the metric function hij (x, y) in scattering coordinates 

(x, y). We shall need to work on the a blown up version of the compactified double 
— 2 

space M in order to do this. 

2 2 

Let M b denote the blown-up compactified manifold consisting of M with the 

boundary (dM) 2 blown up in the sense of ^B] (shown in Figure We define 

p = x'/x" and label the hypersurfaces of M b by lb ('left boundary'), rb ('right 
boundary') and bf ('blown-up face') according as they arise from the faces x' = 0, 

x" = or x' = x" — of M 2 , respectively. Let U be any neighbourhood in M 2 of 
the closure of the diagonal {{z, z) : z G M} U {w £ bf : p(w) — 1, y'(w) = y"(w)}. 
For instance, U could be the region {(z', z") S M : dAi(z', z") < c({z') + {z"))} for 
some small < c <C 1, together with the boundary of this region on the blown up 
face bf . 

Proposition 9.2. Let fl(T], ro) C m\ he defined by 

Q(r],r ) = {x',x" < r^ 1 , d dM (y',y") < v} 
Then for ro sufficiently large and r\ sufficiently small, 

(9.7) the function ^liilfJ ls ln C°°(f7(?7, r ) \ XT). 

Remark 9.3. Notice that U includes a conic neighbourhood of the diagonal near 
infinity, that is, a set of the form 

(9.8) {(r',y',r",y")\ \ r l - \\+ d dM {y' ,y") < e} 
for some e > 0, because of the topology of M b . 

Proof. To prove the claim, we use a symplectic technique; namely, we use the 
function $ = du[z' ^z") as a generating function for a Lagrangian submanifold. 
This function satisfies the eikonal equation 

(9.9) |V Z ,$| 2 = 1, 

and vanishes at the diagonal. Consequently the graph of e?$ will be tangent to the 
flow in SC T*(M) x SC T*(M), lifted to a vector bundle over M b , generated by the en- 
ergy |(|C'lg — !)■ This flow is given by (|9.2|) for the singly-primed variables, while z" 
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diagonal 




Figure 1. The blown up space M b 



and £" are stationary under the flow. Using coordinates (p, y' , x" , y"; v' , //, v" , //'), 
the vector field and the evolution of is given by 



(9.10) 



d 








ds^ 


—x'pv 




—X 

ds 


d ( 'V 

= 






-y 

ds 




■'■'h' 1 , - 


(x') 2 dh^ . , 
f 2 dx>^ 




ds 


-V 

ds 


d , 




x' dh ij , , 




dS^ = 


-x'n' k v' - 


2 d(y') k ^ r 


ds 



ds 



$ = 1. 



Here hij is always a function of the (x',y') variables. Geometrically we consider 
the submanifold L determined by the 'initial condition' that 



(9.11) {x' = x", y> 



y 



// = -m", v' = -v", h^ix'.y')^ + v 2 = 1} C L, 



and L is invariant under the forward flow of (|9.10() (i.e. for positive 'time'). This 
submanifold is a smooth Lagrangian manifold, with boundary, given by the graph 
of o?^, at least near the diagonal. (It is perhaps counterintuitive that L is smooth 
near the diagonal, since <& is singular there. The singularity of $ is reflected in 

2 

the fact that L projects diffeomorphically to M b near the diagonal, but not at the 
diagonal. The inverse image of a point on the diagonal is a (n— l)-sphere consisting 
of all possible unit directions at that point. Smoothness of L is reflected in the fact 
that d<& does not blow up at the diagonal, rather it becomes 'multivalued'.) 

Smoothness of L persists as long as the vector field of which it is the flowout 
remains smooth and nonvanishing. We want to show that the function $/(r' + r") 
extends to a smooth function on (including its boundary at infinity). To do 
this, we shall multiply the vector field l|9.10|l by a factor, which does not change the 
flowout. The new vector field will be smooth and nonvanishing up to the boundary 
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in suitable coordinates, and this will imply that L is a smooth submanifold in these 
coordinates. 

By symmetry of $ under interchange of z' and z", it is sufficient to work in the 
region {r" < 2r'} = {p < 2}. We shall refer to £l(rj, ro) n {p < 2} as the 'region of 
interest'. In the region of interest it is equivalent to require that k = x'<$> = $/r', 
as opposed to $/ (r' +r"), is smooth on Ot, . Its evolution is given by 

d , d d , , , , 

— K = X — < P + <p — X = X — X KV . 

ds ds ds 

As a first step we multiply the flow (|9.1U|I by (a/)" 1 , obtaining the new flow 
d_ 

ds' p -~ PV ± , _ _ , , 1 dhP , , 

(9.12) ±tf)* = h«4 d s^—^-2 8 { y')^ 
ds' J d , 

where ds' := x'ds. This flow is defined for all 'time', but the vector field vanishes 
where p = 0, p = 0, so we need to desingularize about this set (i.e. blow it up). 
Let us introduce the coordinates Mi = p'.J 'p, K = (re — l)/p, and consider the flow 
just on the energy surface {{v') 2 + h 1 ^ p^p'j = 1}. (Introduction of M is equivalent 

to working on the bundle SC T*(M) x SC T*(M) lifted to ~m\ and then blown up at 
{p = 0, p' = 0}). In terms of these quantities, we get, after dropping the primes on 
v and y and changing to a new parameter given by ds" = p ds' 

d 

(9.13) ^ /,<•'. U, ds " 2d V k 
ds" y 3 _d_ _ WMjMj 

±u = P h^M i M j + ^^M i M j '^T- 
ds" 2 dx 

This remains nonsingular as long as \M\ remains bounded and 1 + v ^ 0. For the 
remainder of this paragraph, let us consider what happens at the face bf = \x" = 0}, 
to which (|9.13() is tangent. At bf, the flow is identical to the flow for the perfectly 
conic metric dr 2 +r 2 h(0). Since d s » (h lJ M t Mj) = and initially, -1 < \M\ h < 1, 
M certainly remains bounded. As for v, we have v' > 0, and under the flow H9.13|l . 
d s n p = — v. Thus if p becomes smaller than 1 along a trajectory, then at some 
point v > and therefore v > thereafter along the trajectory. Hence at p = 0, 
v > and there is no singularity in the equation for d s »K. Hence the flow at bf is 
well behaved in the region of interest, and it is not hard to show that it exits the 
region of interest in finite 'time' (measured by the parameter s"). 

By the theory of ODEs depending on a smooth parameter, for sufficiently small 
initial values of x" and for initial values given by l|9.11|l . \M\ remains bounded and 
1 + v strictly positive in the region of interest, and the region of interest is exited 
in finite 'time'. We see then that for small x" , the submanifold L is smooth in the 
region of interest, and K, and hence also k, is a smooth function on L. 

Next we claim that in 0(77, ro)n{p < 2}, for 77 sufficiently small and ro sufficiently 
large, (y',y",p,x") furnish coordinates on L. At L n {x" — 0}, $ is given by the 
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explicit conic formula Ij9.3|l when y' and y" are close: 

$ = ((r') 2 + (r") 2 - 2r'r"cosd dM (y', y ")) 1/2 = r'(l + p 2 + 2pcosd 9M (y' ,y")) 1/2 . 
This is smooth away from 

(9.14) p=l, y'=y". 

Therefore the graph of g?<&, restricted to x" — 0, has (y 1 , y", p) as smooth coordinates 
on any open set excluding i|9.14fl . and hence, their differentials are linearly indepen- 
dent on L n {x" — 0}. Also, the differential of x" is nonzero when restricted to L at 
x" = 0. This follows from the fact that the differential is zero at the initial hyper- 
surface H9.11|) . and that the Lie derivative of dx" with respect to H9.13JI vanishes. 
Therefore, L is projectable — in other words, dy' ,dy" ,dp, dx" have linearly inde- 
pendent differentials when restricted to L at bf, on any open set excluding l|9.14|l . 
and therefore form coordinates on L — uniformly up to the corner {x" = p = 0}. 
By continuity, this remains true on a neighbourhood of L n {x" = 0}, which covers 
an f2(?7,ro) for sufficiently small i] and large r$. Therefore, k is a smooth function 
of (y',y",x",p) on f2(ry, ro) \ U for sufficiently small ry and large rg. Hence for ro 
sufficiently large and rj sufficiently small, there is exactly one geodesic between any 
two points (r',y"),(r",y") with r',r" > r and deuiy' ,y") < V tna t lies wholly 
inside f2(5?7, ro/2). Lemma |9 . 1 1 shows that every globally minimizing geodesic lies 
in this region. Therefore, all geodesies lying wholly within fl(5rj,ro/2) are globally 
minimizing geodesies. This proves that $ = r'n is the genuine distance function on 
M. This completes the proof of the proposition. □ 

We are now able to approximate the asymptotically conic metric g?m rather 
precisely by the perfectly conic analogue (iconic = iconic, h- 

Proposition 9.4. We have 

(9.15) d M (z', z") = d conic (z', z") + e(z', z"), 

where d con i c is the distance function for the perfectly conic metric dr 2 +r 2 hij (0)dy l dyj 
(which is then given by l|9.3|l ). and the error e(z',z") is C°° on ^(77, ro) \ U. 

2 

Proof. We claim that n and K C onic agree at each boundary hypersurface of M b . In 
fact, at lb = {p = 0} they are both equal to 1, while at bf, they solve the same 
ODE with the same initial conditions. Equation (|9.15|l follows immediately. □ 

We can now prove Lemma 16.11 We first prove the claim for perfectly conic 
metrics outside of U. 

Lemma 9.5. Let U be the region U := {(z',z") G M : d M {z',z") < ({z') + 
(z"))/100}. Then we have 

|Vy'Vp"Vv,.«doanic(«',« /, )l < C m ^ m n{ Z ')- m \z")- m " 

for all m', m" > and z' , z" S fl(r/, r )\U . 

Proof. We of course use (|9.3|l . Without loss of generality we may take (z') > (z"); 
since we are outside U, we thus have clonic comparable to (z 1 ). We begin by 
computing derivatives of e? 2 on j c . Observe that cosddM{y\y") is a smooth function 
of dgM{y' ,y") 2 ■ This function is in turn smooth for y',y" G BdM{yo,'n) for V 
sufficiently small (see e.g. J2), so one can then verify the bounds 

IV? V^d 2 oni c(^")l < C m ,, m „^- m \z")^ m " 
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when ml > and m" > 1 (because the (r') 2 term vanishes); for pure z' derivatives 
we have the slightly different formula 

\V^dl ovic { Z ',z")\<C m .{z'f- m '. 

Expanding the latter inequality using the Leibnitz rule, one can then prove 

|V"d conic (z',z")| ^Cm^z') 1 -™' 

inductively for all m! > 0, and then by expanding the former inequality using the 
Leibnitz rule one can also prove 

\VZ'V%"d conic (z',z")\ < C^^iz')-™'^") 1 -™" 

inductively for all m' > and m" > 1. The claim follows. □ 

Next, we verify symbol estimates on a compact portion of the manifold. 

Lemma 9.6. Let zq G M, and let rj > be sufficiently small depending on zq. 
Then we have 

\V% >z ,,d M (z',z")\ < C ZQ , v , m d M {z\z") x - m 
for all m > and z' , z" € B (20,77). 

Proof. We work in normal co-ordinates z 3 around z', so (z') = and J2j W ( z ") \ 1S 
comparable to dM{z', z"). Write r :— d^i(z', z")/ij, and rescale the co-ordinates by 
r (cf. Lemma RT2|) . giving rise to a new manifold M T with metric g 3 T k (z) = gi k (Tz). 
In the ball < C^Ij these metrics g T vary smoothly as r varies over the 

compact set [0, 1], and the distance functions dM T (z', z") also vary smoothly in the 
region where ^ \z^{z')\ < Cr,, £\ \z?(z")\ < C V , and \z\z') - z^(z")\ > cr,. 
In particular we have estimates of the form 

\V%, 2 „d MT (z',z")\ 

uniformly in r in this region. The claim then follows by undoing the scaling. □ 

Now we combine these two estimates to obtain Lemma 16. II 

Proof of Lemma WT\ We let U be the region {(z',z") £ f2(r7,r ) : dj\/(z',z") < 
c(r' + r")} for some small c > 0, together with the boundary of this region on bf . 
Outside of U, we have du(z' ',z") comparable to r' + r", and the symbol bound 
on d]u follows from l|9.15|l . since the error e being smooth on the compactified 
manifold automatically obeys symbol estimates (indeed it obeys estimates which 
are one order of r 1 + r" better than required), while the symbol estimates c? CO nic 
follow from Lemma f9. 51 Inside the region U, we can then rescale as in Lemma 
to rescale z' and z" to a fixed compact set, and then apply Lemma EJlfl if we choose 
c sufficiently small. Note that the claim concerning the nature of the singularity 
of dM at the diagonal z' = z" just follows from the fact that d 2 M is smooth and 
vanishes to second order at the diagonal (see e.g. QU). □ 

It remains to prove Proposition 16.31 We first show this in the case of perfectly 
conic manifolds (in which the non-angular error terms do not appear). 
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Lemma 9.7. For the perfectly conic distance d CO mc{z' , z"), we have the estimate 18 

(9.16) ^- 2 d conic (At), z"(t)) > + J^f) , 

for ddM (y' , y") less than the injectivity radius of (dM,h(0)), where z'(t) moves 
along a geodesic (for the conic metric) with initial condition {z',v'), and z"(t) 
moves along a geodesic with initial condition (z",v"). Here v' ± resp. v'j_ denotes 
the component of v' resp. v" perpendicular to both d r and 4-"f, the tangent vector 
to the geodesic 7 from z' to z" . 

Note that, in the notation of Proposition l6~31 \v' ± \ < |w ang |, \v'[_ \ < \v" ng \, hence 
this result implies Proposition 16.31 in the perfectly conic case. 



Proof. Let us decompose the tangent vectors v' = v' ± + v' pax , v' = v'\_ + Up ar , where 
tip ar is in the span 19 of 7' and d r , and v' ± is perpendicular to the span of 7' and d r , 
and similarly for v" . Clearly the left-hand side of Ij9.16|l is a bilinear form in these 
velocity variables. 

First suppose that v'^ = v'\_ = 0. The geodesic 7(3) lies in a flat plane which 
is totally geodesic. Hence, in this case all geodesies 7(f) from z'(t) to z"(t) lie in 
a flat plane, and the truth of (|9.16|l follows from its truth in Euclidean R 2 , which 
was shown in Section [21 

Next suppose that v' pal = v p ' ar = 0. In this case, we use the formula (|9.3[) and 
compute explicitly, exploiting the fact that in this case 

±d M {z\z")\ t=Q = j t r'(t)\t=o = j/'{t)V=* = j t d 9 M(y',y'%= = 0. 

For a perfectly conic metric we have from (19. 21 that 

^(r') 2 (i)=2|//| 2 , |./ WH/i f/ r ', 

hence d 2 / dt 2 d con i c (z' , z") is given by 
(9.17) 

d- (bo (^ ((/)2) + ^ ((r " )2) - 2(( ^ r> " + { ^ r ' v) C0SddMiy ' y " ] 

+2r'r"^ C o S d 9M (y',y"yj 

Note that (cos ddM (y' , y")) = (1 — dgM(y', y") 2 /2 + 0(dg M j) is a smooth function 
of (y',y"). A computation in normal coordinates around (say) y' shows that 



dt 2 



(co8d eM (y',y")) > -c(^ + lj ^) 2 d d M(y',y"). 



^We caution the reader that there are two 'time' variables in play here; the variable t used 
to parameterize the moving particles z'(t) and z"(t), and also the implicit variable s used to 
parameterize the geodesic connecting z'(t) and z"(t). 

19 If 2' and 2" lie on the same radial arc, then 7' and dr are parallel, but the argument still 
works (e.g. by adjoining an arbitrary vector to the span of d r to define the v' paI component) in 
this case. 
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Since ddM(y',y") is bounded above, we can estimate H9.17(l from below by 

(9 .s,_4_U r+l ,/ f) - 2 (^ + ^)-c r v»(M + ^ ) A 

d conic {z\z") \ r' r" r' r" J 

We rewrite this as 

1 ( 2( (r>-r»)\f'\ 2 . {r"-r')\^ r>f rVf A 

' 4onic0',2")\ ^' r " J 

Since |r' — r"| < d con i c (z' , z"), the first term satisfies the estimate (|9.16|) . For the 
second term, if r' < 2r" < 4r' then it automatically satisfies the estimate (|9.16H . 
while if the ratio of r' and r" is large then d con i c (z', z") is comparable to the larger 
of r' and r", showing that in this case also the second term satisfies the estimate 

Finally consider the case when one vector, say v', is angular, v par = 0, and the 
other is parallel, v'[_ = 0. In this case, we write 

iconic — V( v / )V //)^conic = ^( v ' ; o)^ con i c ^^(v / ) 0)^(0,u // )^conic "i" V(o,t>")^ con i c ' 

and note that the geodesic 7t from z'(0) to lies in the flat plane containing 7 
for all t. Hence v' = v' ± is perpendicular to j t for all t, and we obtain 

V v >d(z',z"(t))=0 foralH. 

Differentiating with respect to t then shows that the cross term in vanishes. 
The estimate H9.16H then follows, since the terms V? v , \d C omc and V? 0it) „)d C onic 
have already been treated. □ 

Lemma 9.8. Suppose that {(r',y'), (r" ,y") are in a small neighbourhood U of the 

diagonal in M b , say in the set (|9.8(l for e > sufficiently small. Then the estimate 
holds. 

Proof. We begin by scaling the problem so that the initial point z'(0) = (r'(0), y'(0)) 
is at r = 1. Let t _1 = r'(0). We scale by factor t; that is, we consider the metric 
g T = dr 2 + r 2 h(r/r, y) and replace the geodesic (r(t), y{t)) by (rr(t), y(t)) which is 
a (reparametrized) geodesic with respect to g T ■ An equivalent formulation of the 
lemma is to prove that 

£_ 

dt 2 ' 



T ti, onic ( z '(t),z"(t)) > -c(K| 2 + K| 2 



-c- 



-{\v'\ 



,."|2 



V 



for all geodesies joining points (r',y'), (r",y") satisfying r' = 1, \r" — 1| < e, 
d dM (y',y") < r). 

To do this, we return to the second variation formula, and use the fact that 

d 2 r d{z'{0),z"{0)) 



dt2 -d M (z'{t),z"{t)) > - J (R(T, J)T, J)ds 

where J is the Jacobi field corresponding to the family of unit speed geodesies 74 (s) 
between z'(t) and z"(t) and T — d^/ds. Let the boundary value of J at z' be 
decomposed as above in to J = J' ± + Jp ar , and define J" and J p ' ar similarly. Let 
Ji be the Jacobi field with boundary values J' ± , J" and let J2 be the Jacobi field 
with boundary values Jp ar , </p' ar - For the perfectly conic metric at r = 0, J2 is in 
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the plane spanned by ^7 and d r . Hence in the general case J2(t) is in the plane 
spanned by ^7 and d r up to 0(t). This means that for any V, 

(R(T,J 2 )T,V) and (R(T, V)T, J 2 ) = 0(t\J 2 \\V\). 

Hence, using the fact that d(z' (0) , z" (0)) is bounded on our rescaled region of 
interest, 

(9.20) ^ d{z '(t), z "(t)) > -0(^1+ r 2 p2\t). 

For r £ [0, 1], we can estimate the length of a Jacobi field corresponding to v = 
(v',v")by 

(9.21) \\J\\l<C(\vf + \v"\ 2 ) 

with a uniform constant C using a comparison estimate as in the proof of Theorem 
4.5.1 of Applying (|9.21|l separately to J\ and J 2 , we find that (|9.2U|) implies 
the desired estimate. □ 

We can finally give the proof of Proposition 16.31 which is the final geometric 
lemma required for Theorem 1 1.51 

Proof of Provosition 1 6'. 31 In a conic neighbourhood U of the diagonal this is pre- 
cisely the content of Lemma 19.81 Outside U, we use Proposition 19.41 and split 
g?m(<z' , z") into the conic distance function d con i C (z' , z") plus a symbolic error e. 
The smoothness of e on Q(r], ro)\U gives an error term of the form 

M 2 W'\ 2 

(9 - 22) J?f + lyf 

As for the conic metric, we need to compare second derivatives of the distance 
function along geodesies for <?, to second derivatives of the distance function along 
geodesies for (/conic- The geodesies differ by 0(t 2 /r 2 ) as one can see from the 
geodesic equation 

df 2 ~ > k dt dt 

since Christoffel symbols for g differ from Christoffel symbols for g con ic by 0(r~ 2 ) 
(with respect to a unit length frame). Hence this discrepancy also gives rise to an 
error of the form l|9.22() . Finally the conic calculation gives rise to error terms of 
the form l|9.16|l . This completes the proof. □ 

10. Appendix II: Proof of local smoothing estimates 

In this section we prove the local smoothing estimates (|l.llf) - (|1.13|) and Lemma 
14.31 The proof of these types of estimates in manifolds usually proceeds by the 
positive commutator method using pseudo-differential operators which are adapted 
to the geometry of geodesic flow. We shall continue to follow this method, however 
instead of using standard pseudo-differential operators, we shall use the calculus 
of scattering pseudo- differential operators on M, described by Melrose ^Bj. This 
calculus is based on the quantization of certain types of scattering symbol classes 
S™' (M), which differ from standard symbols in a number of ways (in particular, 
near dM, every derivative in z of the symbol gains a power of r = X/x), and are 
defined as follows. 
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Definition 10.1. Let m, I be real numbers and < p < 1. A smooth functions 
a : SC T* M — ► C is said to be in the scattering symbol class S™' (M) provided that 
one has the usual Kohn-Nirenberg symbol estimates 

iv«v^( 2 ,c)i <c aAK a + \c\r- lsl 

for z ranging in any compact subset K of M, and all a, j3 > 0, as well as the 
scattering region estimates 

|^V^v£a(z,y,^M)| < C a ,^, s (l + \v\ + i^j)^— y-|*| aj i-«-/»C-y+l*l) 

whenever we are in the scattering region (0, eo) x dM, where the dual variables v, p 
are as in the previous appendix. Equivalently, using radial and angular derivatives 
instead of x and y derivatives, we have 

d?f d2^a(l/r,y,u^) < C a , Pn , s (1 + \v\ + ^Q^-l^-'-^l+^+H). 



This defines seminorms || • || a)/ g )7l <5 on S 1 ™^ (M) in the usual way. Every symbol a in 
S™' (M) can be quantized to give a scattering pseudo- differential operator Op(a) in 
the class \E'™'' ;p (M); the exact means of quantization is not particularly important 
for our purposes, but we can for instance use the Kohn-Nirenberg quantization (on 
any n-dimensional asymptotically conic manifold) 

Op(a)u(z') := (2tt)-" / x(2 , ,z")e 4< - cxp r' 1(2: " ) ^ ) a(z , ,C)w(z")^"dC, 



to define the kernel of Op(a) near dM; here, \ is a cutoff near the diagonal, chosen 
so that the inverse of the exponential map is well defined on supp \- 

Example 10.2. If to, I e R, then an operator of the form (I+H ) m / 2 {z)~ l or (z)~'(l+ 
H) m / 2 lies in \I>™' i;p for any < p < 1. Heuristically, the quantization of a(x, y, /Lt, v) 
can be thought of as a(l/r, y, \d r , iy). 

Remark 10.3. Scattering symbols S™' (M) differ from their more standard coun- 
terparts S™q (M) in that there is a decay of r~ l in the symbol, that differentiation 
in the spatial directions gains a power of r, and that differentiation in the frequency 
directions loses a power of r p . The case p = is the most classical case and is used 
to prove 1)1. 11[) - p. 13(1 : however the proof of Lemma \4 . 31 rea uires the use of more 
exotic symbols with < p < 1 , so that differentiating in the frequency variables v, ji 
costs us a small power of r; this is due to a certain cutoff in p which is necessary 
in our argument. 

We now review the calculus for these scattering pseudo-differential operators, 
obtained 20 in 16 . If A is an operator in *™'' ;P (M), then its symbol a(A) <E 
S™p l (M) is well defined modulo a lower order symbol in S™~ 1 ' l+1 ~ P (M), and its 
adjoint A* is also in this class with symbol 

(10.1) <r(A*) =<7pJ + 0(S"7 U+1 ~ P (M)) 

where we use 0(S™' p l (M)) to denote an error in the class S?J(M). If A e ^ l ' p (M) 
and B € y™'> l ''<P(M) then 

(10.2) AB e *™ +m '' w ' ;p (M); i[A,B] € $w+m'-i,z+z'+i-^(]^ 

2( ^Strictly speaking, these results were only obtained in 1161 for p = but the case for < p < 1 
follows from the same argument. 
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and in fact wc have the more precise formulae 

v(AB) = a(A)a(B) + 0(S^+ m \M)) 

a(i[A,B}) = {a{A),a(B)} + 0(S™+ m '- 2Ml ' +2 - 2p (M)) 

where {/, g} is the usual Poisson bracket on the cotangent bundle T*M. Recall that 
the Poisson bracket {f,g} may also be written as Xf{g) where Xf is the Hamil- 
ton vector field associated to /. This identification is crucial in the commutator 
arguments which follow. 

We introduce the weighted Sobolev spaces H m ' l (M) as 

H m ' l (M) := {u : (z) l u E H m (M)}. 

From itTOI) . ifTOj) and L 2 -boundedness of \l/^ 0;p (M) it is easy to verify that op- 
erators in ^™> l '' p (M) map H m '' l '(M) to H m ' ~ m ' 1 ' +l (M) . 

Now let u(t) evolve via the flow (|1.4|) . Sinced e~ ltH is unitary and commutes 
with (1 + H) s / 2 , we have the trivial estimate 

(10-4) / Ht)\\ 2 H ., 0(M) dt= \\u(Q)\\ 2 Hsm 

Jo 

for any s€R. We now obtain a local smoothing estimate and a Morawetz estimate 
for this flow; note that we are now putting the angz weights into Sobolev norms 
rather than writing them explicitly: 

Lemma 10.4. For any s£R and e > 0, we have 

(10-5) ( |K*)||^,_ 1/2 _ e(M) dt<C S;E ||u(0)||^_ 1/2(M) 

** 

and 

(10-6) J \\ X ?u(t)f He ^ 1/2(M) dt<C s \\u(0)\\ 2 HS _ 1/2{M) 

where x is a smooth cutoff to the scattering region M\Kq. 

Observe that this lemma, combined with the above machinery, immediately gives 

Proof. We first observe that it suffices to prove these estimates for a single value 
of s, since the general case then follows by applying a power of (1 + H) (which 
commutes with e~ ltH ); note that any error terms arising from commuting (1 + H) 
with, for instance, x^ 7 , can be dealt with by l|10.2|l and (|10.4|) . 

We begin by proving the weaker estimate l|1.10fl , which in this language becomes 

(10-7) / IM*)llko ( M) dt<C StV \\u{Q)f He . 1/2{M) 

j o 

for any compactly supported non-negative bump function ip. As before it suffices 
to prove this for a single value of s, say s = 0. This estimate is essentially in 
UJ) |1U| and we only give a sketch here. Using the non-trapping hypothesis of 
M, we may construct a symbol a(z,() £ £]~o' (M) such that {a(H),a} > for 
all ICIg ^ 1 ( sa y)i with the uniform bound {a(H),a} > c > when z is in the 
support of (p. Note that the Poisson bracket {a(H),a} is nothing more than the 
derivative -^a(z(s), ((s)) along the geodesic flow (|9.2|l . which escapes to infinity 
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by hypothesis. See ^D] for the details of this standard construction. If we 
let A — Op(a) G '°''°(M) be the quantization of this operator, then it maps 
H- 1 / 2 (M) to H 1 / 2 (M), and thus 

|(Au(t),«(t)) M | < C\\u(t)\\ 2 H _ 1/2(M) = C||u(0)||^_ 1/a(M) 

for t = 0, 1. Thus by (|2.1() we have 

(i[H,A]u(t),u(t)) M dt < C\\u(0)f H _ 1/2(M) . 

But from the positivity of {a{H), a} we may write 

{a(H),a} =c^+ \e\ 2 + 0(S^\m)) 

where e S 5 1 ' (M) and c > is a constant. Quantizing this using p0.3(l . we 
obtain 

i[H, A] = ap + E*E + O(^^}' -°(M)) 
for some E <G S^'q (M), where 0(^™' l ' ,p [M)) denotes a scattering pseudo-differential 
operator in the class ;p (M). The term E*E gives a positive error which can 
then be discarded, while the contribution of the O( v E'~ 1:0;0 (M)) error is bounded 
by (|10.4|) . and the claim follows. 

Now we prove (|10.6(l . Here the most convenient value of s is s — 1. Let 
a(x,y,^, v) := tp v, where ip{x,y) is a smooth cutoff to the region \x\ < e, and 
< e <C eo is a small number to be chosen later; compare with Example 12.41 Let 
A be the quantization of a. Since a S S 1 ' (M), we have A 6 ^ r ^ 0;0 (M), and we use 
(|2.1(l as before to obtain 

(i[H,A]u(t),u(t)) M < C\\u(0)\\ 2 H1/2{M) . 
From (|9.2|) we see that 



,dh jk 
dx 1 

The second term is dominated by the first in the region < e if e is sufficiently 
small, 
write 



ih 2 dh? k 

region \x\ 

small, while the third term is compactly supported and in S 1 ' (M). We can thus 



{a(H), a} = c^xh^fi^k + \e\ 2 + p 2 0(Sl'°(M)) + 0(S$(M)) 

for some e 6 SI'q(M) and some compactly supported function <p(z). Quantizing 
this we obtain 
i[H,A] = c^y^)*,^^^^*^^ 

The second term is again positive and can be discarded, while the third term can 
be controlled using l|l(J.7fl and the fourth term by I|1U.4|) . The claim <|1(J.6[) follows 
(the contribution of the region |x| > e being controlled by (|10.7I) ). 

Finally we prove (|1U.5|) . Again we set s = 1, and now use the symbol a := 
—ip 2 x 2£ v where ip is as before. Since a S ^^(M), the quantization A is in 
*^ 0;0 (M), so by (flTT) as before we have 

" (i[H,A}u(t),u(t))M<C\\u(0)f H1/2(My 
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From (H3 we have 

{a(H), a} = 2ex 1+2e v 2 ^ 2 - x 2e {a(H), ifv}. 
From the proof of l|l(J.6fl we see that 

{<j(H)^ 2 v} = x^^^OiS^iM)) + ^ 2 O(S 2 : ° (M)) + O(Sl:° (M)) 
for some compactly supported tp, and thus 

i[H,A] = Ce(x 1/2+e d r TP)*{x 1/2+£ d r iP) 

+ f^jWpf)^ + pO(^°(M))p + 0(^°(M)). 

The second error term is controlled by (|10.6|l . the third by (|1U.7|) . and the fourth 
by l|10.4f> . This proves l|10.5|l in the region |x| < e; the region |x| > e is of course 
controlled by ifTITTjl . □ 

To prove Lemma 14.31 it turns out not to be practicable to apply the positive 
commutator method directly, mainly because the z integrals in Lemma 14.31 are 
not themselves positive. Instead, we use the positive commutator method to first 
estimate an auxiliary positive quantity, which we then use to control the expressions 
in Lemma l4.3l More specifically, we prove the following variant of H10.6J) which only 
uses 'half an angular derivative instead of a full angular derivative. 

Lemma 10.5 (Half-angular Morawetz estimate). Let < p < 1, and let <fi ■ R. — ► R 
be a smooth non- decreasing function such that <f> = on (— oo, 1] and <f> = 1 on 
[2, +oo). Let b(x,y, p,,v) be the S\ ,X J 2 {M) symbol 

(10.8) b := l^m 2 + k| 2 )0(^)- 1/2 (0(^)IH 1/2 |H 1/2 + iMl) 

where |/i| 2 = ImI^x) := ^ k ( x ^y)t JL jf J 'k- Let B € ^^^^{M) be a quantization of 
this operator. Then for any H s+1 / 2 solution to 11.4JI and any sdR, we have 

(10.9) f \\Bu(t)\\ 2 Hs(M) dt < C\\u(0)\\ 2 HS+1/2{M) . 
Here the constants depend on p and s. 

Remark 10.6. Ignoring all the cutoffs, the important part of b is the l/xj 1 / 2 ^! 1 / 2 
term. Heuristically, B can be thought of as the operator r -1 / 2 !)? 7 ! 1 / 2 ^,! 1 / 2 , with the 
various cutoffs involving <f> needed to avoid a singularity arising from any degeneracy 
of y or V, localizing physical space to the scattering region r > 1/eo, and frequency 
space to the region \ J7\ > |V|r~ p > r~ p . This cutoff is not dangerous for us as in 
the region where |y| < |V|r~ p , one can effectively deduce Lemma f4.3l from (|10.5f) 
(this is why we require p > 0). The reader may wish to ignore the presence of 
the 4> cutoffs for a first reading. The estimate H10.9J) can then be thought of, very 
non-rigorously, as 

f \\xM 1/2 \f\ 1/2 u(t)\\ 2 H^ HM} dt < C|K0)|| 2 ff3+1/2(M) 

J 

for some appropriate cutoff %; compare this with (|10.(i|) . The crucial point here is 
that we do not lose an epsilon of decay as we would from l|10.5fl . 
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Proof. As in the proof of Lemma, ri0.4l it suffices to verify this when s = 0. We again 
use the positive commutator method. Morally speaking, the commutant to use is 
|y 7 |sgn(i9 r ); however this operator is too singular, and so we must use a modified 
version of this commutant. 
Let a denote the function 

a{x,y,n,v) := 2 (e o /x)0V + | M | 2 ) (V(i)0 2 (j^j)(-HsgH + Cv 

where C is a large constant to be chosen later. Notice that a vanishes in the near 
region K . This can easily be verified to be a scattering symbol of order S 1 '(M), 
because of all the cutoffs. Let A <G ^lf' p (M) be the quantization of a. Applying 
(|2.1|) as in Lemma Tl 0.41 we have 

" (i[H,A}u(t),u(t))M<C\\u(0)\\ 2 Hl/HM y 
We shall now factorize i [H, A] in the form 

(10.10) i[H,A] = B*B + C*d + 0(^ 1+5;p (M)) + 0{^°' P (M)) 

finite 

where C t € $k 1/2 '' p (M) and S > 0; this will prove the claim since the error terms 
can be treated by l|10.5|l , l|10.4fl ■ From (|10.3|) it suffices to obtain a representation 



(10.11) Wfl),o} = \b\ 2 + V M 2 + 0{S 2 ^ p +s {M)) + o(sl : ° p (M)). 



Observe that {a(H), a} vanishes when x > eo, and in the region eg > x > eo/4 



it is a symbol of order Sy. P (M) as mentioned earlier, and hence in this region 
is also a symbol of order S 2 ' l p (M) for any I. Thus this factorization is easy to 
obtain in the near region, and we can focus on the scattering region x < eo/2. In 
particular we can now ignore the cutoff (f>(eo/x). A similar argument allows us to 
work in the region where \/i\ 2 + \u\ 2 > 2 (as {a(H), a} is certainly in S 1 ' p (M) when 
|/i| 2 + \v\ 2 < 4), allowing us to ignore the </>(|/-t| 2 + \ v \ 2 ) cutoff. We thus need to 
consider 

MH), ^ 2 (i^)0 2 (M) HA1 | sgm ,) + Cv) 

in the region x < eo/2, |^| 2 + \v\ 2 > 2. 

Since the two 4> 2 terms are each 1 on the support on the derivative of the other, 
we may express this as 

(^){a{H)-Wgnv}+C{a{H),v} 



(10.12) HX " ,M| 

-|/i| S gw{ ( T(iJ),0 2 (M) } _| A1 | sglw{(7(jff);0 2 ( H )}. 

We have {a(H),-\fi\sgnu} = x(\n\\v\ + 0(5 2 j(M))) and {a(H),v} = x{\^\ 2 + 
0(S 2 ;*(M))). Hence for C > 1 the sum of these terms is larger than \b\ 2 , where b 
is as in (110. 8(1 . To analyse the third term we compute 

MiJ),0 2 (g )} = 2!B #'(M)(| M | + ^)+ 0(S 2 p 2 (M)). 
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This term is dominated by the second term, modulo acceptable errors, for suffi- 
ciently large C. As for the fourth term, we compute 
(10.13) 

MWD, 2 (^)} - -2x| M | # '(-^L)((i±|M + M. + 0(( ^(M))) . 

x p \v\ xP\v\ V x p x p v z ' H I 

These terms are all supported where \p\ is comparable to x p v. Since there is an 
overall factor of |/x| 2 on the right hand side of (|1U.13|) . we gain a factor of x 2p , 
showing that (|1U.13[1 is a symbol of order (2, 1 + p) which is an error term in 

Assume that the function <f> is such that terms such as \/(f> and ^Jlj? ', etc, are 
smooth; for example, we make take 

<p(t) = 1 - 0(3 - 1) and <f>(t) = e- (t - irl for t G [1, 1.2]. 

Then {a(H),a} can be written, modulo acceptable errors, as \b\ 2 plus a sum of 
squares \°i\ 2 ■ This completes the proof of IjlU.lUp . which establishes the lemma. 

□ 

Proof of Lemma \4 ■ °J[ We first prove 1)4.5(1 . Observe if the derivative was re- 
placed by an angular derivative then this claim would follow from i|1.13|l or i|10.6fl . 

so it suffices to show that 



sup 

o wew 



a J w (z)fMt,z)d r u{t,z) 
dg(z) 



M 



dt 



<C\\u(0)f H1/2(M) 



where a 3 w satisfies similar symbol estimates to a{f. Choose a < p < 1 (for instance, 
we may take p := 1/2). Observe that we may write 



ai ] {z)y j u{t,z)d r u{t,z) f 

dg(z)= A w u(t, z)u(t, z) dg(z) 

<M r J 

where A w € ^ 2 ^ P (M) has symbol 

a(A w )(x,y,u, p) = Ca? w {x,y)xp j v + 0(S['°) 

for some constant C (indeed we may improve the error term substantially). It is 
here that we require the functions a J w to be symbols; otherwise A w would not be 
a pseudodifferential operator. Observe that in the region \p\ 2 + \v\ 2 < 4, the main 



term is in . From the support hypothesis of a 3 we thus have 



o-{A w )(x,y,v,p) = C(j){e /x)(j){\p\- + \v\ )a : > w {x,y)xp Jj v + 0{S 1 \ p ). 

Recalling the operator B G *sc 1/2;p (M) with symbol b G S\^J 2 (M) fr om the pre- 
vious lemma, we can thus factorize 

a(A w )(x,y, V ,fx)=O(S°l p (M))\b\ 2 

+ C^eo/xMtf + |H 2 )(1 - 4>{^))<{x, y)xp,v + 0(Slf p (M)). 

The second term is supported on the region where \p\ < 2\v\x p and can thus be 
easily seen to lie in S{ p +P (M). Thus by (|1U.3|) we have 

A w = B*0(^' p (M))B + 0(^ 1+p:p ) + 0{^lf' p ); 
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since V™* 1 * maps H m / 2 ^ l / 2 {M) to H- m / 2t / 2 (M), 

| J A w u{t,z)u{t,z)dg{z)\ <q| J B^)|| L2 (M)+C|| U (i)|| ffl ._i+ £(M) +C'|| U (<)|| H i (M) 

uniformly in w. The claim (|4.5|l then follows from Lemma TlU. 51 (|l(J.5fl and l|l(J.4|l . 

Now we prove (|4.6|l . As before we have 

f a{ u (z)f i u(t,z)u(t,z) f _ 

/ dg(z) = / A w u(t, z)u(t, z)dg(z) 

Jm r J 

where A w € *^ 1;p (M) has symbol 

cr(A u ,)(a;,y,^ //) = Ca : > w {x,y)x^ j + 0(5°'°), 

which is similar to the previous but without the factor of v. As before, we localize 
the main term to the region |/x| 2 + |^| 2 > 1. We can then smoothly split the main 
term into the region where \v\ > and where > In the first region, we 
have (1 + a(H)) 1 / 2 comparable to and so by repeating the previous argument 
we eventually obtain 

A w = B*O('9~ c 1 ' '' p (M))B + 0(* S V + " ;P ) + O(*°< 0;p ), 

and one then argues as before. In the second region, we have (1 + <j{H)) 1 / 2 com- 
parable to |/x|, so we can write directly 

a{A w ){x,y,p,n) = O{S^ p )x\^\ 2 + 0(5°'°), 

and hence we have a decomposition 

A w = yo(* s - 1 '°^(M))y + o(^) 

uniformly in w. The claim now follows using i|10.6|) instead of Lemma ll0.5l □ 
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